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1.  INTRODUCTION 


1.1  Objective 

Tfce  overall  objective  of  this  project  is  to  develop  a  computer 
code  which  simulates  the  performance  characteristics  of  a  hydrazine- 
driven,  turbo-alternator  power  system  with  operating  time  (propel¬ 
lant  storage  capacity),  type  of  propellant,  and  power  output  as  in¬ 
dependent  variables.  Within  the  scope  of  this  part  of  the  project, 
operating  time  (two  hours)  and  type  of  propellant  (hydrazine)  re¬ 
mained  fixed.  Power  output  was  the  system  independent  variable  and 

the  effects  of  varying  power  output  on  the  entire  system  were  in¬ 
vestigated. 

1.2  Explanation  of  Problem  and  Approach 

Ihe  turbo-alternator  power  system  consists  of  a  small,  multi¬ 
stage,  gas  turbine  which  drives  a  high  rotational  speed  alternator. 
Schematic  and  pictorial  diagrams  are  shown  in  Figs.  1  and  2,  re¬ 
spectively.  Preliminary  calculations  and  available  data  indicated 
that  a  turbine  powered  with  the  products  of  catalytic  decomposition 
of  liquid  hydrazine  has  certain  advantages  over  other  propellant 
systems.  Liquid  hydrazine,  in  the  presence  of  a  suitable  catalyst, 
decomposes  to  form  hydrogen,  nitrogen,  and  ammonia  (the  amount  of 
ammonia  produced  is  a  function  of  temperature  and  decomposition 
chamber  residence  time).  Hydrazine  was  chosen  as  the  propellant  to 
be  modeled,  while,  bearing  in  mind,  other  mono-propellant  and  bi¬ 
propellant  combinations  could  be  added  to  a  future  simulation  program. 

The  system  simulation  must  be  detailed  enough  to  insure  that 
reasonably  accurate  operating  characteristics  can  be  predicted.  This 
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necessitates  use  of  accurate  propellant  properties  and  detailed 
analysis  of  the  decomposition  chamber  and  gas  turbine.  The  proposed 
program  began  at  a  simple  stage  with  greatly  simplified  equations 
with  many  limiting  assumptions  and  proceeded  in  a  stepwise  fashion 
to  a  gradually  more  complex  system  of  equations  with  fewer  limiting 
assumptions.  Throughout  this  program,  an  operating  time  of  approxi¬ 
mately  two  hours  was  assumed  reasonable.  For  this  operating  time 
a  propellant  storage  capacity  of  approximately  14.75  m3  is  required. 
TVo  hours  is  also  a  safe  upper  .limit  in  light  of  present  state-of- 
the  art  technology  of  hydrazine  decomposition  chambers  and  small, 
high-speed  gas  turbines. 


2.  GENERALIZED  SYSTEM  SIMULATION  PROGRAM 

2.1  Turbo- Alternator  System  Simulation 

A  system  is  a  collection  of  related  components,  and  in  the 
turbo-alternator  system  the  components  considered  were  the  de¬ 
composition  chamber,  the  turbine  stages,  and  the  exhaust  system- 
Ihe  term  "system  simulation"  may  be  defined  as  observing  a  syn¬ 
thetic  system  that  closely  imitates  the  performance  of  a  real 
system.  As  a  system  simulate,  this  program  predicts  the  steady- 
state  operating  quantities  within  the  system  (pressures,  tempera¬ 
tures,  energy-  and  fluid -flow  rates)  at  the  condition  where  all 
energy  and  material  balances,  all  equations  of  state  of  working  sub¬ 
stances,  and  all  performance  characterise . .  of  individual  components 
and  controls  are  satisfied.  In  the  case  of  the  turbo-alternator,  no 
real  system  yet  exists.  The  equations  used  in  the  simulation  are  all 
either  fundamental  equations  developed  for  any  similar  turbine  power 
system  or  empirical  equations  fit  to  data  determined  for  similar 
system  components  (i.e.  ,  decomposition  chamber).  Ihe  equations  are 
nearly  all  non-linear  in  nature,  and,  as  a  result,  some  special  simul¬ 
taneous  method  of  solution  was  required.  The  method  chosen  was 

"Newton -Raphson"  and  a  particular  generalized  system  simulation  program 
was  used.  See  Ref.  [23].* 

2.2  Newton-Raphson  Technique  for  more  than  One  Equation 

Trie  Newton -Raphson  technique  is  an  iterative  method  for  solving 
n  simultaneous  non-linear  equations  for  n  unknowns.  Ihe  general  steps, 

^Numbers  in  brackets  refer  to  entries  in  REFERENCES. 
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the  background  for  which  is  found  in  Ref.  [24],  are: 

1.  Move  an  the  terms  to  one  side  of  each  equation, 

yi(xl*x2*-*-»Xn>  =  0 

yn(xi*x2»‘*‘»xn)  =  0 

2.  Assume  trial  values  of  the  unknowns  x  - 

It*  “2t*‘"  * 

V  and  substitute  into  the  equations.  Ain***  certainly 

the  values  of  through  yn  win  „ot  be  2ero>  so  caU  these 

values  R  through  R  . 

1  n 

3.  The  corrections  to  be  made  ar*>  fntmr?  ,  . 

raae  are  found  by  solving  the  following 

set  of  simultaneous  linear  equations  for  the  Ax's. 

R  s  A  ^i 

1  X1  +  x2  +  •••  +  Axn 

\  =  r2-  Ax,  t  Ax  +  +  . 

n  1  ^7  2  *  •  •  +  Axn 

£  n 


where  ^  .  W.  ^  ..  , 


’2  ,new 


)»  etc. 


*’  ”  ^  =hangeS  “  iU  °f  the  *’•  «•  sufficiently  small,  the 

computation  can  be  terminated.  If  not.  ^  tales  cn  the  value 

l,new*  2t  becomes  x2 .new *  stc-»  and  the  process  returns 
to  step  2. 

In  review,  the  liewton-Raphson  method  converts  the  solution  of  a  sat  of 
n  Simultaneous  nonlinear  equations  to  an  iterative  process  each  iteratic 
Of  which  requires  the  solution  of  a  set  of  n  simultaneous  linear  equatioi 

2.3  Generalized  Computer  Program 

H.e  generalized  program  used  was  taken  from  [2],  It  „as  chosen  be. 
cause  of  this  author's  familiarity  with  both  the  trthod  and  the  program 
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aid  the  program’s  adaptability  to  this  system.  From  the  standpoint  of 
computer  calculations,  the  operations  are  as  follows: 

1.  Introduction  of  trial  values  of  variables, 

2.  Computation  of  R’s, 

3.  Calculating  the  partial  derivatives, 

4.  Solving  the  simultaneous  linear  equations , 

5.  Checking  the  changes  in  the  x's  against  some  convergence 
criteria,  and 

6.  If  convergence  is  not  satisfied,  compute  improved  values  of  the 
x's  and  return  to  operation  2. 

Operation  1  is  the  only  unique  one  when  moving  from  one  system  to 
another,  although  to  execute  operations  2  and  3  a  subroutine  must  Le  avail¬ 
able  which  can  be  called  to  supply  the  equations  unique  to  the  system  at 
hand. 

A  listing  of  the  generalized  program  is  provided  in  the  appendix. 

The  complete  program  consists  of  (1)  the  main  program,  (2)  the  equation 
subroutine  EQNS,  (3)  the  subroutine  for  extracting  the  partial  deriva¬ 
tives  PARDIF ,  (4)  the  simultaneous  linear  equation  solver  GAUSSY,  and 
(5)  the  data  cards.  The  main  program  first  reads  in  the  data  cards  and 
immediately  prints  out  the  values  that  appear  on  these  data  cards.  Next 
the  program  initializes  the  iteration  counter,  because  a  feature  of  the 
program  is  that  if  a  specified  number  of  iterations  are  performed,  the 
program  terminates  even  though  convergence  may  not  yet  be  achieved.  The 
program  then  calls  EQNS  subroutine  to  compute  the  values  of  the  R's. 
Following  this,  the  main  program  calls  PARDIF  which  in  turn  calls  EQNS 
to  numerically  determine  the  partial  derivatives.  The  values  of  the 
derivatives  may  be  of  some  interest,  particularly  for  purposes  of 
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analysis  if  the  program  fails.  Most  of  the  partial  derivatives  will 
be  zero,  especially  in  large  systems,  so  only  the  non-zero  derivatives 
are  printe  out.  The  partial  derivatives  and  the  values  of  R  are 
transferred  to  the  linear  equation  solver  GAUSSY.  The  printout  now 
provides  the  results  of  this  iteration.  Finally,  the  number  of  iter¬ 
ations  an  the  closing  tolerance  on  the  change  of  all  variables  is 

checked,  and  the  program  either  terminates  or  returns  for  another 
iteration . 
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3.  TURBO-ALTERNATOR  SYSTEM  AND  SYSTEM  COMPONENTS 


3.1  Decomposition  Chamber 

The  decomposition  chamber  presents  a  great  challenge  in  terms 
of  simulation.  The  processes  Involved  in  the  catalytic  decomposition 
of  hydrazine  are  complicated  and  extremely  difficult  to  accurately 
model.  Upon  reviewing  jet  propulsion  laboratory  abstracts,  and  the 
STAR  (Scientific  and  Technical  Aerospace  Reports)  index,  ten  different 
papers  and  reports  were  found.  The  most  useful  among  these  was  a  United 
Aircraft  Corporation  report  [22]  which  included  a  computer  program  for 
one-  and  two-dimensional  simulations  of  catalyst  chambers.  The  program 
itself  was  not  used  in  this  simulation;  however,  several  useful  empiri¬ 
cal  correlations  were  obtained  from  the  United  Aircraft  studies.  The 
two  equations  which  are  represented  graphically  in  Fig.  3  and  4  are 
taken  directly  from  typical  catalyst  chamber  data.  The  equations  repre¬ 
sent  functional  relationships  between  chamber  pressure,  chamber  tempera¬ 
ture,  mass  flow  of  hydrazine,  average  catalyst  particle  size,  length  and 
diameter  of  cylindrical  chamber,  and  fractional  dissociation  of  ammonia 
gas.  These  equations  are  applicable  primarily  in  cases  where  most  of 
the  hydrazine  decomposition  occurs  in  the  first  few  millimeters  of  the 
reactor;  this  rapid  hydrazine  decomposition  rate  is  associated  with  re¬ 
actors  packed  with  particles  25  mesh  or  smaller  for  approximately  5  mm 
of  bed  length.  For  these  cases  the  correlations  work  well  for  axial 
distances  greater  than  25  mm  and  for  values  of  pressure  between  69,000 
and  690,000  N/m2  (10  to  1000  psia),  a  mass  flux  between  7,05  and  70.5 
kg/m2-sec  (1,44  to  14.4  lbm/ft2-sec) ,  and  equivalent  spherical  radius 
between  0.3  and  3  mm.  For  a  reactor  packed  with  small  (<  25  mesh) 
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Temperature  in  Interstitial  Phase  T(z),  °R 


Figure  3  Empirical  Correlation  for  Temperature  in  a 
Hydrazine  Decomposition  Chamber 


10 


Figure  4  Empirical  Correlation  for  fractional  A-monia  Dissociation 
in  a  Hydrazine  Decomposition  Chamber 


11 


particles  for  the  first  few  millimeters  and  larger  particles  thereafter, 
the  particle  radius  refers  to  the  larger  particles. 

In  order  to  use  these  correlations  in  this  program,  it  was  neces¬ 
sary  to  determine  the  approximate  chamber  dimensions  for  such  an  appli¬ 
cation  and  an  estimate  of  catalyst  particle  size  was  also  needed.  In 
gas  generation  applications,  the  peak  performance  is  achieved  with 
values  of  fractional  dissociation  of  ammonia  of  0.60  to  0.70.  By  using 
this  criteria  and  running  a  simple  program  with  various  values  of 
chamber  dimensions,  a  set  of  reasonable  dimensions  were  determined. 
Throughout  the  entire  project  the  following  equipment  parameters  of 
the  cylindrical  decomposition  chamber  were  fixed:  z  =  chamber 
length  =  0.203  m,  D  =  chamber  diameter  =  0.254  m,  and  A  =  catalyst 
particle  radius  =  3.048  mm. 

3.2  Turbine 

Since  the  hardware  for  this  system  is  presently  non-existent,  a 
great  deal  of  work  was  involved  in  the  turbine  modeling.  It  was  not 
only  necessary  to  develop  the  simulation  equations  but  also  to  do  a  pre¬ 
liminary  design  so  that  the  model  could  be  considered  representative  of 
a  typical  system.  Initial  design  criteria  involved  the  use  of  dimension¬ 
less  flow  coefficients  and  stage  loading  coefficients  to  determine  the 
approximate  stage  configurations.  Because  of  small  inlet  areas  and  small 
blade  heights,  it  was  necessary  to  use  partial  admission,  impulse  staging 
for  the  first  stage.  The  second  and  third  stages,  which  were  de¬ 
termined  necessary  to  meet  the  power  requirements,  seemed  best  suited  to 
50  percent  reaction  staging  because  of  the  greater  efficiencies  attained. 
The  blade  velocities  were  assumed  not  to  exceed  520  m/sec  (1700  ft/sec), 
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this  being  based  upon  presently  available  materials  which  can  withstand 
the  stresses  developed  at  such  speeds.  Design  curves  and  criteria  were 
taken  frcm  Horlock  [9],  Shepherd  [21],  Balje  [3],  and  Duisinberre  [4]. 
With  the  aid  of  these  curves  and  the  equations  upon  which  they  were 
based,  physical  dimensions,  angles,  and  dimensionless  parameters  were 
determined  for  the  design  load  of  2  mW.  Many  of  these  parameters  vary, 
of  course,  at  off -design  conditions.  All  of  these  off -design  parameters 
were  maintained  as  variables  and  may  be  determined  for  any  conditions  by 
the  computer  program. 

The  modeling  of  the  turbine  consists  of  three  distinct  section::: 
stages  1,  2,  and  3.  Each  section  requires  five  equipment  parameters: 

(1)  ALF  =  angle  of  gas  exit  velocity  vector  from  nozzle,  (2)  BET  = 
angle  of  relative  gas  exit  velocity  vector  from  rotor,  (3)  U  =  blade 
speed,  (4)  A  =  nozzle  throat  area,  and  (5)  HOB  =  ratio  of  blade  height 
to  length.  All  gas  properties  and  velocities  are  calculated  in  each 
stage.  The  specific  heats  and  the  specific  heat  ratios  of  the  gas  are 
calculated  at  a  mean  temperature  and  pressure  in  each  stage.  Loss  char¬ 
acteristics  based  upon  blade  angles  and  velocity  vector  angles  were 
taken  from  Horlock  [9].  More  specifically,  losses  are  calculated  using 
the  Soderburg  loss  coefficients  for  both  the  stator  and  the  rotor.  Tne 
total-to-total  efficiency  of  each  stage  which  uses  total  (stagnation) 
values  of  pressure  and  temperature  as  references  is  a  function  of  these 
loss  coefficients,  the  axial  velocity  of  the  gas,  the  blade  speed,  and 
the  various  inlet  and  outlet  angles.  The  power  produced  at  each  stage 
is  taken  from  the  definition  of  expansion  efficiency  where  the  expans:on 
efficiency  is  the  ratio  of  actual  work  to  the  work  developed  in  an  isen- 
tropic  expansion.  The  actual  work  is  obtained  through  the  loading 
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coefficient  relationships,  while  the  isentropio  work  is  obtained  through 
simple  isentropic  pressure  ratio  relationships. 

The  total  power  developed  by  the  turbine  is  simply  the  sum  of  the 
powers  developed  by  each  of  the  three  stages.  A  total  of  62  equation, 
was  required  to  satisfactorily  model  this  turbine.  The  stage  equations 

"*  MSily  bK*en  <iC*m  lnto  groups  of  equations  representing 

3  Sl"Ele  SU8e  “*•  «  is  -  -T  operation  then  to  modify  the  turbine 

design  by  adding  or  subtracting  stages  or  changing  any  of  the  five  equip- 
aent  parameters  associated  with  a  stage. 

3.3  Exhaust  System 

The  exhaust  system  consists  simply  of  an  exhaust  duct  terminating 
at  a  converging  nozzle.  The  nozzle  is  to  be  choked  at  all  operating 
conditions.  For  the  purposes  of  this  simulation,  the  exhaust  system 
was  assumed  to  obey  isentropic  flow  characteristics.  The  duct  cross- 
sectional  area  was  set  at  0.00811  „2;  however,  the  throat  area  was  al¬ 
lowed  to  change  as  a  variable.  The  reasons  behind  allowing  the  throat 
area  to  change  lie  in  system  design;setting  the  throat  area  presents  a 
constraint  upon  the  system  which  must  he  satisfied,  and  at  this  stage 
final  system  parameters  have  not  been  determined.  It  was  easier  to  keep 
throat  area  as  a  variable  and  note  bow  it  changed  at  different  load  con¬ 
ditions.  It  was  necessary,  of  course,  to  set  the  throat  velocity  in 
order  to  insure  choked  conditions.  In  the  final  system  simulation, 
throat  area  would  likely  be  an  equipnent  parameter  which  would  not 
change.  This  modification  will  easily  re  implemented  in  the  present 


program. 
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3.4  Variables  Used 

System  or  Equipment  Variables 

=  effective  rotor  exit  area  of  a  stage;  i*2,4,5  for 

2 

stages  1,  2,  3,  respectively,  m 
C.  =  gas  velocity  at  ith  stage  nozzle  exit;  i  =  1,  3,  5  for 
stages  1,  2,  3,  respectively,  m/sec 
Cj  =  gas  velocity  at  jth  stage  rotor  exit;  j  =  2,  4,  6  for 
stages  1,  2,  3,  respectively,  m/sec 
C^.  =  axial  velocity  in  ith  stage;  i  =  1,  2,  3,  m/sec 
CP^  =  mean  specific  heat  of  gas  through  ith  stage;  i  =  1,  2,  3, 
J/i2-°K 


CP^  =  mean  specific  heat  of  exhaust  gas,  J/kg-°K 

2 

GO  =  mass  flux  through  decomposition  chamber,  kg/sec-m 
KGi  =  mean  specific  heat  ratio  of  gas  through  ith  stage; 
i  *  1,  2,  3 


KG^  =  mean  specific  heat  ratio  of  exhaust  gas 
m  =  mass  flow  rate  of  fuel,  kg/sec 
MWT  =  molecular  weight  of  ouik  gas,  kg/mole 

=  static  pressure  at  a  stage  inlet;  i  =  1,  3,  5  for 
stages  1,  2,  3,  respectively,  bars 
P.  =  static  pressure  at  a  stage  rotor  exit;  j  =  2,  4,  6  for 

J 

stages  1,  2,  3,  respectively,  bars 

P  .  =  total  pressure  at  a  stage  inlet;  i  =  1,  3,  5  for 
oi 

stages  1,  2,  3,  respectively,  bars 
P^7  =  total  pressure  at  third-stage  rotor  exit,  bcxs 
Pp  =  static  pressure  in  exhaust  duct,  bars 
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P,p  =  static  pressure  at  exhaust  nozzle  throat,  bars 
POWi  =  power  output  of  ith  stage,  i  =  1,  2,  3,  m 
Ti  =  static  temperature  at  a  stage  inlet ;  i  =  1,  3,  5  for 
stages  1,  2,  3,  respectively,  °K 
Tj  =  static  temperature  at  a  stage  rotor  exit;  j  =  2,  4,  6 
for  stages  1,  2,  3,  respectively,  °K 
TQi  =  total  temperature  at  a  stage  inlet;  =  1 ,  i ,  5  for 
stages  1,  2,  3,  respectively,  °K 
Tq7  =  total  temperature  at  a  third-stage  rotor  exit,  °K 
TD  =  static  temperature  in  exhaust  duct,  °K 
T,j,  =  static  temperature  at  exhaust  nozzle  throat,  °K 
TRSO  =  turbine  power  output,  mw 

VD  =  gas  velocity  in  exhaust  duct,  m/sec 
VT  =  gas  velocixy  at  exhaust  nozzle  throat,  m/sec 
XH2  =  mo^e  fracti°n  of  hydrogen 
XJJ2  =  mole  fraction  of  nitrogen 
=  mole  fraction  of  ammonia 

cm  ~  angle  of  actual  gas  velocity  vector  at  a  stage  nozzle 
exit;  i  =  3,  5,  7  for  stag'  >  1,  2,  3,  respectively 
0^  =  angle  of  relative  gas  velocity  vector  at  a  stage  rotor 
exit;  i  =  2,  4,  6  for  stages  1,  2,  3,  respectively 
nTTi  =  totai"to~totai  efficiency  for  the  ith  stage;  i  =  1,  2,  3 
Pi  =  gas  density  at  a  stage  inlet;  i  =  1,  3,  5  for  stages 

3 

1,  2,  3,  respectively,  lg/m 

=  gas  density  at  a  stage  rotor  exit;  j  =  2,  4,  6  for 

3 

stages  1,  2,  3,  respectively,  kg/m 
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| 

PD  -  gas  density  in  exhaust  duct ,  kg/m3 
PT  s  gas  density  at  exhaust  nozzle  throat,  kg/m3 
$  =  fractional  dissociatiru  of  ammonia 
5Ni  *  Soderburg  loss  coefficient  for  ith  stage  statorj 
i  =  1,  2,  3 

^Ri  =  Soderburg  loss  coefficient  for  ith  stage  rotor; 
i  =  1,  2,  3 

System  of  Equipment  Parameters 

A  =  average  spherical  catalyst  particle  radius,  m 

\  *  area  °f  ith  stage  nozzle  exit,  i  *  1,  3,  5  for  stages  1, 

2 

2,  3,  respectively,  m 

Aq  =  cross  section  area  of  exhaust  duct,  m2 
At  =  cross  section  area  of  exhaust  nozzle  throat,  m2 
ALTO  =  alternator  power  output ,  mw 

D  =  diameter  of  decomposition  chamber,  m 
HOBi  =  ratio  of  blade  height -to-length  in  ith  stage; 
i  =  1 ,  2 ,  3 

Ui  =  rotov  blade  velocity  in  ith  stage;  i  =  1,  2,  3,  m/sec 
Z  =  length  of  decomposition  chamber,  m 
°i  =  an®^e  °f  actual  gas  velocity  vector  at  a  stage  nozzle 
exit;  i  =  2,  4,  6  for  stages  1,  2,  3,  respectively 
3^  =  angle  of  relative  gas  velocitj  vector  at  a  stage  rotor 
exit;  i  =  3,  5,  7  for  stages  1,  2,  3,  respectively 
nA  =  alternator  efficiency 
hg  =  gearbox  efficiency 
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4.  GAS  PROPERTIES 


4,1  Decomposition  Process 

Ihe  gas  properties  and  their  accurate  determination  is  one  of 
the  nos.-  important  factors  in  this  type  of  simulation  model.  An  ex¬ 
tensive  literature  search  was  carried  out,  and  a  great  deal  of  useful 
information  .ibout  the  gas  properties  was  found  in  Reid  [18],  McBride 

[15],  Ellenwood  [5],  Hilsenrath  [8],  Kubin  [11],  Obert  [16],  and 
Van  Wylen  [25], 

The  actual  hydrazine  decomposition  was  assumed  to  be  a  two-step 
process : 

1.  The  reaction  of  hydrazine  to  ammonia,  hydrogen,  and  nitrogen 

2N2H4  =  2“H3  +  »2  +  N2  (4.1) 

?•  Ihe  dissociation  of  ammonia  into  hydrogen  and  nitrogen 

2NH3  =  N2  +  3H2  (4.2) 

The  preceding  equations  are  the  same  equations  used  by  United  Air¬ 
craft  Corporation  [22]  in  their  analysis  of  the  phenomenon.  One  im¬ 
portant  assumption  was  made  at  this  point:  Cbce  the  gas  leaves  the 
decomposition  chamber,  no  more  ammonia  dissociates  and  the  bulk  gas 
molecular  weight  remains  constant  through  the  remainder  of  the  system. 
Actually,  this  assumption  is  quite  valid  and  is  substantiated  by  ex¬ 
perimental  data.  Prom  Eqs.  (4.1)  and  (4.2),  it  is  a  simple  matter  to 

determine  all  coefficients  in  terms  of  the  fractional  dissociation  of 
ammonia. 

N2H4  =  (1  '  ^  NH3  +  J  (1  -  <P)  H2  +  I  (1  +  3(t>)  H2  (4.3) 
w.  are  <j>  is  the  fractional  dissociation  of  ammonia.  From  Eq.  (4.3),  the 
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mole  fraction*  of  each  gas  can  be  expressed  as 
tional  dissociation  of  ammonia. 

a  function  of  the  frac 

Xjijj  2  (l  ~  <fr)/(2  +  Q) 

3 

(4.4) 

X«  *  (1  +  4>)/( 4  +  241) 

2 

(4.5) 

Xj|  *  (!  +  34>V(4  +  241) 

(4.6) 

where  X  is  the  mole  fraction  of  each  gas  in  the  bulk  gas  mixture.  The 
molecular  weight  of  the  bulk  gas  is  related  to  the  gas  mole  fractions 
and  the  gas  molecular  weights  through  Eq.  (4.7). 


MWT  =  *NH  *17-032  +  Xj,  *28.016  +  X  *2.016  (4.7) 

3  2  °2 

4.2  P-p-T  Relations 

Tbe  previous  equations  indicate  the  assumption  of  ideal  gas  behavior. 
Actually,  the  gas  was  assumed  to  be  semi-perfect  in  behavior;  that  is,  the 
relationship  P  =  pRT  was  applied,  however,  certain  thermodynamic  properties 
were  assumed  to  be  functions  of  temperature  and  pressure.  The  applic¬ 
ability  of  ideal  gas  laws  was  based  upon  calculations  of  deviation  from 
ideal  behavior  using  compressibility  factors  (see  APPENDIX  A).  Although 
these  calculations  indicate  only  deviations  of  1  to  3  percent  from  ideal 
gas  behavior,  the  deviations  for  specific  heats  and  specific  heat  ratios 
of  the  gases  were  much  greater. 

4.3  Specific  Heat  Relationships 

Zero  pressure  specific  heat  data  for  each  of  the  cases  were  found  in 
McBride  [15],  Hilsenrath  [8],  and  Kubin  [11].  Third-order  polynomial  ex¬ 
pressions  for  these  specific  heat  functions  of  temperature  were  found  in 
Oben  [16]  and  7an  Wylen  [25],  ir.ese  equations  were  converted  to  System 
International  units,  and  specific  heat  values  were  generated  for  different 
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temperatures  and  checked  with  the  other  data  sources  mentioned  above. 

The  agreement  between  the  calculated  and  documented  values  was  excel¬ 
lent.  The  one  minor  drawback  to  these  equations ,  however,  was  their 
absence  of  pressure  corrections.  As  is  shown  graphically  in  Figs.  5 
and  6  (Ellenwood  [5]),  there  is  a  dependence  upon  pressure ,  especially 
at  lower  temperatures .  The  effects  of  pressure  upon  specific  heat  as 
shown  in  Figs.  5  and  6  were  determined  using  the  Beattie -Bridgeman 
equations  of  state.  This  method  of  specific  heat  calculation,  although 
quite  accurate  and  fundamental,  is  cumbersome.  Rather  than  trying  to 
use  Beattie-Bridgeman  equations,  the  polynomial  expressions  for  zero 
pressure  specific  heats  were  "adjusted"  to  take  into  account  the  small 

variations  in  the  regions  of  534°K  to  1368°K  (500°F  to  2000°F)  and  0  to 
2 

13,800,000  N/m  (0  to  2000  psia).  The  three  specific  heat  equations  were 
combined  through  the  gas  mole  fractions  to  give  the  specific  heat  of  the 
bulk  gas  at  any  temperature  and  pressure  in  the  above  mentioned  region. 

The  ideal  relationship  between  specific  heat  at  constant  pressure 
and  specific  heat  ratio  was  used  to  calculate  specific  heat  ratio. 
(Equations  showing  these  relationships  will  appear  in  following  chapters.) 
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Figure  5  The  Effect  of  Temperature  on  C  of  Nitrogen  at  Various  Pressures 


Effect  of  Temperature  on  C  of  Hydrogen  at  Various  Pressures 


5.  APPLICATION  OF  FUNDAMENTAL  EQUATIONS 


The  equations  used  to  model  the  system  are  described  below.  The 
second  stage  of  the  turbine  was  not  specifically  included  because  it 
is  the  same  as  the  third  stage.  (The  sources  and  derivations  of  these 
equations  may  be  found  in  APPENDIX  A.) 


5.1  Alternator  Power  Equation 


_  _  ALTO 

TBB0  -  TnT^np- 


where 


TRBO  =  turbine  power  output,  mw 
ALTO  =  required  alternator  output,  mw 


=  alternator  efficiency 


rig  =  speed  reducer  (gearbox)  efficiency 


5.2  Turbine  Power  Equation 

TRBO  =  P0W1  +  P0W2  +  P0W3 


where 


TRBO  =  turbine  power  output,  watts 
P0W1  =  power  produced  in  first  stage,  .nw 
P0W2  =  power  produced  in  second  stage,  mw 
P0W3  =  power  produced  in  third  stage ,  mw 


5.3  Mass  Flux  Equation 


;o  = 


4(m) 


•TD 
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(5.1) 


(5.2) 


(5.3) 


I 


fail  ff  riiBiifi 'll 


where 


(5.4) 


o 

GO  =  mass  flux  through  decomposition  chamber*  kg/sec-m 
m  =  mass  flow  rate  of  fuel,  kg/sec 
D  =  diameter  of  cylindrical  chamber,  m 

5.4  Fractional  Dissociation  of  Ammonia  Equation 

/  68.94757  \  0,22  1  1 

*HH  J’0-17} 


where 

4>  =  fractional  dissociation  of  ammonia 
Z  =  length  of  decomposition  chamber,  m 
A  =  average  spherical  catalyst  particle  radius,  m 
Pq1  =  total  pressure  in  decomposition  chamber,  bars 

5.5  Decomposition  Chamber  Temperature  Equation 

/  P 

T  -  I1"20  lU  '  *  [°-°75  (-68TW57-)]}  *  1535  ) 

Ol  .  „ 


(5.5) 


1.8 


where 


Tq1  =  total  temperature  in  decomposition  chamber,  °K 


5.6  Mole  Fraction  of  Ammonia  Equation 

XNH  =  (1  -  <J>)/(2  +  $) 

3 


(5.6) 


where 


X.TI,  =  mole  fraction  of  ammonia 
bn- 
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5.7  Mole  Fraction  of  Nitrogen  Equation 

=  (1  +  <>)/4  +  20) 

"2 


wnere 


yN  *  *»le  fraction  of  nitrogen 
5.8  Mole  Fraction  of  Hydrogen  Equation 


where 


XH2  =  (1  +  30>/(4  +  20) 


XH2  =  mole  fr>action  of  hydrogen 


5.9  Molecular  Height  Equation 


where 


(5.7) 


(5.8) 


MKT  .  (17.032)  *  ,  (28.016)  ^  .  (2.016)  X  (5.9) 

o  o  n_ 


MWT  =  molecular  weight  of  bulk  gas,  kg/mole 

Equations  (5.10)  through  (5.30)  represent  the  set  of  equations  aesorih- 
ing  the  first  impulse  stage . 

5.10  First-Stage  Inlet  Total  Pressure  Equation 


P,  Cf 

Pol  =  pi  *  -V1  10-5 


(5.10) 


where 


Pi  =  Static  pressure  at  first  stage  nozzle  exit,  bars 
°1  =  gas  density  at  first  stage  nozzle  exit,  kg/n3 

C1  =  g3S  Velocity  at  first  stage  nozzle  exit,  m/sec 


5.15  Velocity  Vector  Relationship  Equation 

tan  (a2)  +  tan  (<*3)  =  tan  (62)  +  tan  (8  )  (5.15) 

where 

a2  =  ^g16  of  actual  gas  velocity  vector  at  first  stage 
nozzle  exit,  radians 

a2  =  angle  of  actual  gas  velocity  vector  at  first  stage 
rotor  exit,  radians 


26 


5.16 


5.17 


where 


5.18 


5.19 


where 


=  angle  of  relative  gap  velocity  vector  at  first  stage 
nozzle  exit,  radians 

(?3  -  angle  of  relative  gas  velocity  vector  at  first  stage 
.  rotor  exit ,  raaians 


Impulse  Stage  Angle  Equation 


(5.16) 


First-Stage  Specific  Work  Equation 


2CAX1  V1  ttan  ^3^  =  nTTl  CP1  Tol 


/p  .[(KG  -D/CKG,  )] 

.“(d  ' 

(5.17) 


CAX1  =  axial  velocity  of  gas  through  first  stage,  m/sec 
U1  =  first  stage  blade  velocity,  m/sec 
r,TTl  =  total'to_totai  efficiency  for  first  stage 
Pq3  =  total  pressure  after  first  stage ,  bars 

First-Stage  Power  Equation 


P0W1  =  »  nTT1  cPl  Tol 


[ .KG  -1)/(KG  )] 

■  -  d 


10 


“e  (5.18) 


First-Stage  Exit  Total  Temperature  Equation 


Tc3  =  Tol 1  1  '  nTTl 


/p  [KG  -1)/(KG  )] 

^d 


(5.13) 


Tq3  =  total  temperature  after  first  stage,  °K 
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t>.20  lirst-Stage  Axial  Velocity  Equation 

CAX1  =  C1  Ccoc  <0‘2)] 


(5.20) 


5.21  Rotor  Loss  Coefficient  Equation 


2ir 


V  1  °-°25 


where 


;  ♦(*£)_ 


i  + 


3.2 


HCB1 


(5.21) 


^R1  =  SoderburS  loss  coefficient  for  first  stage  rotor 


H0B1  =  first  stage  ratio  of  blade  height-to-blade  length 


5.22  Stator  (Nozzle)  Loss  Coefficient  Equation 

2. 


N1 


=  o.025  1  -1  -  ][i  +  J.’2 

\  1.57  J  J|^  H0B1 


(5.22) 


where 


^N1  =  Soderbur2  loss  coefficient  for  first  stage  stator 


5.23  First-Stage  Total-to-Total  Efficiency  Equation 

1 


TT1 


(5.23) 


[tan  (a2)  +  tan  (a  )] 


5.24  First -Stage  Exit  Total  Pressure  Equation 

J2 


P  C 

P03  =  P2  *  15 


(5.24) 


where 


Po3  =  static  pressure  at  first  stage  rotor  exit,  bars 
P2  =  density  at  first  stage  rotor  exit,  kg/m3 
C2  =  velocity  at  first  stage  rotor  exit,  m/sec 
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-{[29,087.17  -  l,Wp V  4.00U6  ,  10-3 

'(^)  '  8-88854  x  0'7(^ 

'  (isr)]  \}] f  mT 


(5.29) 


Equation  (5.29)  is  used  several  times  throughout  the  program  to  cal¬ 
culate  a  mean  specific  heat  over  a  small  range  of  temperature  and  pr 
sure  variation. 

5.30  First-Stage  Mean  Specific  Heat  Patio  Equation 


KG  =  - ±__ 

1  CP1  -  R 


(5.30) 


Equation  (5.30)  is  used  to  calculate  a  mean  value  of  specific  heat  ratio 

based  upon  the  specific  heat  at  constant  pressure  calculated  in  Eq.  (5.29). 

Equations  (5.31)  through  (5.50)  represent  the  set  of  equations  used  to 

"Cdei  the  third  stage.  The  second  stage  is  described  by  a  similar  set  of 

equations.  Both  the  second  and  third  stages  were  assumed  to  be  50  percent 
reaction  stages. 

5.31  Third -St age  Inlet  Total  Pressure  Equation 


where 


P5Cr  , 

Po5  =  P5  * 


1 o5  total  Pressure  after  second  stage,  bars 
P5  =  static  pressure  at  third  stage  nozzle  exit,  bars 
P5  -  gas  density  at  third  stage  nozzle  exit,  kg/m3 
C5  =  gas  velocity  at  third  stage  nozzle  exit,  m/sec 


(5.31) 


=  static 
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5.32  Third-Stage  Inlet  Total  Temperature  Equation 


where 


5.33 


5.34 


where 


5.35 


where 


5.36 


ob 


T  +  — 

5  2CP„ 


(5.32) 


Tq5  =  total  temperature  after  second  stage,  °K 


Tg  =  static  temperature  at  third  stage  nozzle  exit,  °K 
CP3  =  mean  specific  heat  in  third  stage,  J/kg-°K 


Third  Stage  Inlet  Equation  of  State 

■,-5 


P5  =  p5  rt5  10- 


(5.33) 


Third-Stage  Inlet  Continuity  Equation 


m  =  p5  C5  A5 


(5.34) 


A  =  area  of  third  stage  nozzle,  m 


Flow  Coefficient  Equation 
U3 


r  =  tan  (3  )  -  tan  (6  ) 
CAX3  7  6 


(5.35) 


U3  =  third  stage  blade  velocity,  m/sec 
^AX3  =  ax^a^-  velocity  of  gas  in  third  ttage,  m/sec 

=  angle  of  relative  gas  velocity  vector  at  third- 
stage  rotor  exit,  radians 

3g  =  angle  of  relative  gas  velocity  vector  at  third- 
stage  nozzle  exit ,  radians 


50  Percent  Reaction  Stage  Angle  Equation 


(5.36) 
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MMiki 


■  iWijf.fttf.ninr 


where 


*/  =  311616  of  actual  S*5  velocity  vector  at  third 
stage  rotor  exit,  radians 


5.37  Third-Stage  Specific  Work  Equation 


where 


U3  CAX3  [tan  (“6>  +  »■>  («,)!  =  rlTr3  CP3 


C(KG3-1)/(KG3)] 


(5.37) 


a6  =  angle  of  actual  gas  velocity  vector  at  third 
stage  nozzle  exit,  radians 
nTT3  =  total-to-total  efficiency  for  third  stage 
Po7  =  total  Pressure  after  third  stage,  bars 
KG3  =  mean  sPeclflc  heat  ratio  in  third  stage 

5.38  Third-Stage  Power  Equation 


P°“3  '  "  TT3  CPi  »ToS 


1  - 


?o7  \  C(KG3"1)/(KG3)]", 

^ j 


>  -6 


10  ”  (5.38) 


5.39  Third-Stage  Exit  Total  Temperature  Equation 

’P  _\  t(KG3-D/(KG  )]n 


ro?  ^os 


1  - 


o7 

P  , 
o5/ 


where 


To7  =  total  temperature  after  third  stage,  °K 

5.40  Third-Stage  Axial  Velocity  Equations 
CAX3  =  C3  Ccos  (a6)] 


(5.39) 


(5.40) 
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5.41  Rotor  Loss  Coefficient  Equation 

5R3  ‘  °-°“  [l  ♦(%*£)  ](X  *  M) 

where 


(5.41) 


£  =  Soderburg  loss  coefficient  for  third  stage  rotor 

R3 

H0B3  =  third  stage  ratio  of  blade  height-to-blade  length 


5.42  Stator  (Nozzle)  Loss  Coefficient  Equation 

2- 


N3 


0.025  1  .(^1)  (l*M) 


(5.42) 


where 


^N3 

“5 


=  Soderburg  loss  coefficient  for  third  stage  stator 
=  angle  of  actual  gas  velocity  vector  at  second  stage 
rotor  exit,  radians 


5.43  Third-Stage  Total-to-Total  Efficiency  Equation 

1 


TT3 


1  +*< 


’N3  \U3 ) 

+  E 

CAX3 

R3 

U3  cos  (0?) 

r 

"i 

>(-w) 

|  tan  (o6)  +  tan  (o^) 

(5.43) 


5.44  Third-Stage  Exit  Total  Temperature  Equation 

p  i 

Po7  =  P6  +  V-  1°'5 


where 


(5.44) 


P.  =  static  pressure  at  third  stage  rotor  exit,  bars 
o 

3 

P-  =  density  at  third  stage  rotor  exit,  kg/m 

b  ~ 

C  =  velocity  at  third  stage  rotor  exit,  m/sec 
6 
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5.45  Third-Stage  Exit  Total  Tenperature  Equat 


ion 


=  T,  + 


07  "  *6  T  2CPJ 


(5.45) 


where 


T6  =  static  temperature  at  third  stage  rotor  exit* 


°K 


5.46  Third- Stage  Exit  Equation  of  State 


P6  *  h  RT6  I0'5 


(5.46) 


5.47  Third-Stage  Exit  Continuity  Equat 


ion 


■  5  p6  C6  “ 


(5.47) 


where 


A6  =  calculated  effective  area  of  third  stage  rot 


or  exit*  m 


5.48  Third-Stage  Rotor  Exit  Velocity  Equati 

C 


on 


C„  = 


AX3 


6  cos 


T^7 


(5.48) 


non 


5.49  Third-Stage  Mean  Specific  Heat  Equati 

This  equation  is  identical  to  Eq.  (5.29)  except  for  different 
n^ean  temperatures  and  pressures: 


CP3  =  f(T5’T6’P5’V 


(5.49) 


5.50  Third- Stage  Mean  Specific  Heat  Ratio  Equat 

CP„ 


ion 


KG„  = 


3  CP3  -  R 


(5.50) 


Equations  (5.51)  through  (5.61)  represent  the  equations  used  to  model 
the  exhaust  system. 
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fur  Hr  nl nil  ri 


5-5i  W“",t  Buct  Total  Pressure  Equation 


where 


P..V? 


<•  =  p  +  _ D  P  ,_-5 

o/  D  5 —  i0 


PoV  =  total  pressure  after  third  stage,  bars 
P“  =  Statis  F^^ure  i„  exhaust  duct,  bare 
PB  "  gas  density  Fn  exhaust  duct,  kg/m3 
VD  *  ges  velocity  in  exhaust  duct,  m/sec 

5-52  Exhaust  Duct  Total  Temperature  Equation 


(5.51) 


T  =  T  + 
07  j£  r 


where 


To7  =  total  temperature  after  third  sta 
Td  =  static  temperature  in  exhaust  due 
CPEX  =  1116311  specific  heat  of  exhaust  gas 

'.53  Exhaust  Duct  Continuity  Equation 


stage,  °K 
duct,  °X 
S°s,  J/kg-°K 


®  *  Pr,  V_  A 
»  i)  D 


where 


aD  =  «ea  of  exhaust  duct  cross  section,  mJ 

5'5“  ::':;iaus*  “oct  Equation  of  state 
Pn  =  Pn  RT„  10-5 


D  -  PD  10  ~3 

5.55  Exhaust  Nozzle  Total  Pressure 


Equation 


_  pt  v£ 

Po7  PT  +  2 —  10~5 


(5.52) 


(5.53) 


(5.54) 


(5.55) 


where 


rT  "  static  pressure  at  nozzle  throat,  bars 
PT  =  gas  density  at  nozzle  throat,  kg/m3 
Vf  =  S35  velocity  at  nozzle  throat,  m/sec 

5.56  Exhaust  Nozzle  Total  Temperature  Equation 


T  =  T  +  _ L 

09  T  2CP 


where 


tt  =  static  tenperature  at  nozzle  throat, 
5.57  Exhaust  Nozzle  Continuity  Equation 


m  =  pT  VT  At 


where 


At  =  area  of  exhaust  nozzle  throat. 


m 


5.58  Exhaust  Nozzle  Equation  of  State 


PT  =  PT  RTt  10'5 


5.59  Exhaust  Nozzle  Throat  Velocity  Equation 

VT 

where 


RT„ 


KGX  =  mean  specific  heat  ratio  of  exhaust  gas 

5.60  Exhaust  Gas  Mean . Specific  Heat  Equation 

Equation  (5.60)  is  the  same  as  Eq.  (5.29)  and  Eq.  (5.49) 

it  is  evaluated  at  a  different  mean  temperature 

CP,  =  f (t  T  P  P  ) 

X  x 1 T* rD,r T 


36 


(5.56) 


(5.57) 


(5.58) 


(5.59) 


except 

(5.6C) 


-•  ■  . . 


MMMiaajkdklajwli 


(5.61) 


37 


6.  DISCUSSION  OF  POSSIBLE  FUTURE  REFINEMENTS  OR  ADDITIONS 


6.1  Decomposition  Chamber 

The  decomposition  chamber  representation  in  the  present  program  is 
one -dimensional  and  cannot  account  for  different  catalyst  bed  or  inlet 
port  configurations .  A  future  change  to  the  program  could  be  to  add 
the  previously  mentioned  Uhited  Aircraft  Corporation  catalytic  chamber 
two-dimensional  simulation  [22].  This  addition  would  definitely  be  a 
better  chamber  representation  and  also  would  be  an  excellent  first 
step  toward  transient  system  simulation.  Ultimately,  in  order  to  model 
a  transient  system,  control  devices  must  be  added.  The  control  system 
models  must  be  complete  enough  to  insure  that  accurate  system  response 
could  be  determined. 

The  effect,  of  inlet  fuel  temperature  on  system  performance  was  not 
considered  in  this  program.  This  effect,  although  small  for  small  in¬ 
let  differences,  should  be  considered  for  the  sake  of  completeness.  A 
complete  discussion  of  the  relationships  between  inlet  fuel  temperatures 
and  chamber  performance  appears  in  Wrobel  [26]. 

6.2  Gas  Properties 

Although  the  gases  behave  nearly  ideal  in  the  operating  ranges  of 
this  program,  some  error  is  introduced  through  the  use  of  ideal  gas 
equations.  These  small  errors  could  be  essentially  eliminated  through 
the  use  of  Beattie -Bridge man  or  Virial  coefficient  equations  of  state. 
The  implementation  of  these  equations  could  become  quite  cumbersome, 
but  the  increased  accuracy  of  the  system  model  may  justify  their  use. 

A  discussion  of  the  Beattie -Bridge man  equation  and  the  necessary 


coefficients  may  be  found  in  Van  Wylen  125],  Qbert  [16],  and  Ellenwood 
[5]. 


5.3  Separate  Turbine  Simulation 

Because  of  the  extensive  nature  of  the  turbine  simulation,  it  may 
become  necessary  to  reduce  the  actual  number  of  modeling  equations. 

This  could  be  accomplished  by  completing  a  separate  turbine  simulation 
program  and  then  using  the  results  of  such  a  program  to  obtain  empirical 
performance  curves .  These  curves  could  provide  accurate  simulation 
equations,  but  far  fewer  equations  would  be  required. 
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APPENDIX  A 

The  following  explanations  are  supplied  for  certain  of  the  unique 
equations  found  in  the  simulation  program.  Many  equations  which  are 
fundamental  (continuity,  stagnation  pressure,  etc.)  will  not  be  discussed. 

A.l  Equations  (5.4)  and  (5.5) 

Equations  (5.4)  and  (5.5)  are  empirical  correlations  based  upon  a 
significant  amount  of  experimental  data.  These  equations  are  repre¬ 
sented  graphically  in  Figs.  5  and  6  taken  from  Smith  [223. 

A. 2  Equations  (5.7),  (5.8),  and  (5.9) 

These  gas  mole  fraction  'nuations  were  determined  according  to 
the  decomposition  reaction  as  presented  in  Eq.  (4.2)  of  this  report 
and  in  Smith  [22]. 

A. 3  Equation  (5.13) 

Equation  (5.13)  is  essentially  a  perfect  gas  correlation.  The  use 
of  perfect  gas  laws  is  substantiated  by  the  following  data  and  calculations. 
Compressibility  factor  charts  were  used  to  determine  the  deviation  of  the 
gases  (under  typical  conditions)  from  ideal  behavior. 

Data 

:  Tc  =  59 . 7°R  and  Pc  =  188  psia 

N2:  Tc  =  227. 0°R  and  Pc  =  493  psia 

NH  :  T  =  731. 4°R  and  P  =  1657  psia 

j  c  c 

Compressibility  factor  =  Z  =  Fv/RT;  as  Z  -*■  1  the  gas  approaches  ideal  be¬ 
havior.  Then  (assuming  typical  extremes  of  T  =  1500°R  and  P  ■■  1500  psia) 
from  the  generalized  chart  (Fig.  A.l,  p.  621  [25]),  for  H0,  Z  =  1.05; 
for  ,  Z  =  1.01;  and  for  NH^,  Z  =  0.99.  These  compressibility  factors 
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show  that  deviation  from  1.0  (ideal  gas)  is  very  small  and  as  temperatures 
become  higher  and  pressures  become  lower,  the  chart  shows  even  closer  to 
ideal  behavior  is  observed. 

A. 4  Equation  (5.15) 

This  equation  relates  the  actual  and  relative  inlet  and  exit  velocity 
vector  angles  for  the  first  two  impulse  stages.  This  relationship  holds 
for  any  degree  of  reaction  and  was  taken  from  Horlock  [9j. 

A. 5  Equation  (5.16) 

Taken  from  Horlock  [9],  this  equation  results  from  the  definition 
of  an  impulse  (zero  reaction)  stage. 

A. 6  Equation  (5.17) 

Equation  (5.17)  is  a  combination  of  the  definition  of  stage  loading 
coefficient  in  Horlock  [9]  and  the  efficiency  of  a  turbine.  Tne  ideal 
work  is  taken  from  the  definition  of  an  isentropic  expansion  and  the 
actual  work  is  derived  from  the  stage  loading  coefficient:. 

A. 7  Equations  (5.18)  and  (5.19) 

Both  Eqs.  (5.18)  and  (5.19)  -ere  developed  directly  from  the  definition 
of  expansion  efficiency  and  specific  work  in  Shepherd  [21], 

A. 8  Equation  (5.20) 

This  equation  is  determined  directly  from  the  velocity  vector  diagram 
of  a  typical  turbine  stage. 

A. 9  Equations  (5.21)  and  (5.22) 

These  two  equations  define  the  loss  coeff'cients  for  a  turbine  stage. 
These  loss  coefficients  developed  by  Soderburg  were  taken  from  Horlock  [9]. 


A. 10  Equation  (5.23) 

Equation  (5.23)  is  the  total-to-total  efficiency  of  a  turbine  stage 
in  terns  of  Soderburg  loss  coefficients,  characteristic  stage  velocities, 
and  angles.  This  equation  was  taken  from  Horlock  [9]. 

A. 11  Equation  (5.29) 

This  equation  for  the  specific  heat  of  the  gas  at  any  mean  tempera¬ 
ture  and  pressure  was  developed  through  the  use  of  polynomial  expressions 
for  the  specific  heats  of  each  of  the  component  gases  given  in  Obert  [16]. 
Each  of  the  component  gas,  molar  specific  heats  is  multiplied  by  the  gas 
mole  fraction  to  obtain  the  proper  "bulk  gas"  specific  heat.  There  is  a 
correction  factor  included  in  each  gas  polynomial  expression  to  account 
for  pressure  effects  on  specific  heat.  These  correction  factors  were  de¬ 
termined  using  Figs.  5  and  6  as  the  basis  for  nitrogen  and  hydrogen,  re¬ 
spectively.  Without  any  data  available,  the  ammonia  factor  was  assumed 
to  be  the  same  as  the  nitrogen  factor. 

A. 12  Equation  (5.35) 

Equation  (5.35)  is  the  definition  of  a  flow  coefficient  for  any 
stage  given  in  Horlock  [9].  This  equation  could  be  used  interchangeably 
with  Eq.  (5.15)  for  a  turbine  stage  in  this  program.  Both  provide  the 
necessary  geometric  constraints  on  the  velocity  vector  relationships. 

A. 13  Equation  (5.36) 

This  equation  is  the  geometric  relationship  which  must  be  satisfied 
in  a  50  percent  reaction  turbine  stage  (taken  from  Horlock  [9]). 
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APPENDIX  B 


*•**•  TURBO-AL TERN# 1 OR  SIMULATION  PROGRAM  WITH  A  3  STAGE  TURBINE  *>**♦* 

c 

c 

c 

C  MAIN  PROGRAM  FOR  GENERALI  'l-.it  S YSTE“  SIMULATION 
C 

t  SUBROUTINES  CONSIST  OF  11)  EQUATIONS  -  12)  PARTIAL  DERIVATIVES  — 

C  13)  GAUSSIAN  ELIMINATION  FOR  SIMULTANEOUS  SOLUTION  OF  LINEAR  EPS 
C 

C  GLOSSARY  OF  TERMS  USED  IN  MAIN  PROGRAM 
C  ITER  a  NUMBER  OF  ITERATIONS 

C  ITMA*  a  MAXIMUM  NUMBER  OF  ITERATIONS  TO  BE  PERMITTED 

C  NVAR  *  NUMBER  OF  UNKNOWNS  =  NUMBER  OF  EQUATIONS 

C  POU.J)  =  PARTIAL  DERIVATIVE  OF  FUNCTION  I  WITH  RESP  TO  VARIABLE  J 
C  Rl  )  *  RESIOUAt.  OF  EQUATION 

C  TLRNCE  *  MAXIMUM  FRACTION  OF  VALUE  OF  VARIABLES  PERMITTED  BEFORE 
C  ITERATION  COMPLETE.  THUS  TLRNCE  *  O.Ol  REQUIRES  CHANGE  OF 

C  ALL  VARIABLES  TO  BE  LESS  THAN  1  PERCENT  FUR  CONVERGENCE 

C  V(  )  *  VALUE  OF  THE  VARIABLE 

C  VCORR (  )  *  CORRECTION  IN  THE  VARIABLE  DURING  THIS  ITERATION 
C  OESf  )  a  DESIGNATION  OF  VARIABLE 

c 

C  INPUT  FORMAT 

C  FIRST-  CARO  -  NVAR  IN  13  FIELD  RIGHT  JUSTIFIED 

C  SECOND  CARO  -  TLRNCE  IN  F10.4  FIELD  FOLLOWED  BY  IT MAX  IN  iO-SPACES 

C  THIRD  AND  FOLLOWING  CARDS  ARE  THE  TRIAL  VALUES  OF  THE  VARIABLFS 

C  IN  8F1G.3  FIELDS  UNTIL  COMPLETE 

C  LAST  SET  OF  CARDS  CONTAINS  VARIABLE  DESIGNATIONS  IN  SEQUENCE  IN 
C  20A4  FIELDS  UNTIl  COMPLETE 
C 

c 

C  SYSTEM  OR  EQUIPMENT  PARAMETERS 

C 

C 

c  A=AVERAGE  SPHERICAL  CATALYST  PARTICLE  RADIUS  (METERS) 

C  AlaAREA  OF  1ST  STAGE  NOZZLE  EXIT  (SQ.M) 

C  '  A3= AREA  OF  2ND  STAGE  NOZZLE  EXIT  (SQ.M) 

C  A5=AREA  OF  3RD  STAGE  NOZZLE  EXIT  (SQ.M) 

C  ADUCT=CROSS  SECTIONAL  AREA  OF  EXHAUST  DUCT  I  SQ.M) 

C  ALF2aANGLE  OF  ACTUAL  GAS  VELOCITY  VECTOR  AT  THE  1ST  STAC-E  NOZZLE 
C  EXIT  (RADIANS) 

C  ALF4= ANGLE  OF  ACTUAL  C-AS  VELOCITY  VECTOR  AT  THE  2ND  STAGF  NOZZLE 
C  EXIT  (RADIANS) 

C  ALF6aANGLE  OF  ACTUAL  GAS  VELOCITY  VECTOR  AT  THE  3RD  STAGE  NOZZLE 
C  EXIT  (RADIANS) 

C  ALTEsALT ERNATOR  EFFICIENCY 

C  ALTOsALTERNATOR  POWER  OUTPUT  (MW) 

C  AT=CROSS  SECTIONAL  AREA  OF  EXHAUST  NOZZLE  THROAT  (SQ.M) 

C  BET3aANGLF  OF  RELATIVE  GAS  VELOCITY  VECTOR  AT  THE  1ST  STAGE  ROTOR 
C  EXIT  (RADIANS) 

C  8ET5»iNGLE  OF  RELATIVE  GAS  VELOCITY  VECTOR  AT  THE  2ND  STAGE  ROTOR 
C  EXIT  (RADIANS) 

C  BET7=ANGL f  OF  RELATIVE  GAS  VELOCITY  VECTOR  AT  THE  3RD  STAGE  ROTOR 

C  EXIT  (RADIANS) 

C  D*OI AMETER  OF  DECOMPOSITION  CHAMBER  !M) 

C  GRRF=GEARflOZ  EFFICIENCY 

C  HOBlaRATIO  OF  BLADE  HEIGHT  TO  LENGTH  IN  1ST  STAGE 
HOB2»RATIO  OF  PLACE  HEIGHT  TO  LENGTH  IN  2NP  STAGE 
H0B3»RATI0  OF  BLADE  HEIGHT  TO  LENGTH  IN  3RD  STAGE 
UlaROTDR  BLADE  VELOCITY  IN  1ST  STAGE  (M/SEC) 
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U2-ROT0R  BLADE  VELOCITY  IN  2N0  STAGE  (M/SEC) 
U3«ROTOfl  BLADE  VELOCITY  IN  3RD  STAGE  IN/ SEC) 
ZslFNGTH  OF  DECOMPOSITION  CHAMBER  CM) 


DECOMPOSITION  CHAMBER  VARIABLES 


VC  1 )*TUR8INE  OUTPUT  (MW)CBRITISH  UNI TS=HOPSEPOWER) 

VC  2)  sHASS  FLOW  PATE  OF  FUEL  CKG/SEOC  BRITISH  UNI  TS*LBM/SEC) 

VC  3)*H4SS  FLUX  THROUGH  DECOMPOSITION  CHAMBER  CKG/SEC-SO.M) 
(BRITISH  UNIT  S=LBV/ SEC -SO.FT ) 

VC  4) rFRACT ION  OF  AMMONIA  DISSOCIATED 
VC  5)=MOLE  FRACTION  OF  NH3  GAS 
VC  6 ) =MOl F  FRACTION  OF  N2  GAS 
VC  7) *MOl F  FRACTION  OF  H2  GAS 

VC  8)=MOFCULAR  WEIGHT  OF  GAS  (KG/HOLE  )  C BRI TISH  UNIT S=LBM/MOLE ) 


TURBINE  VARIABLES 


FIRST  STAGE 


VC  9  )  =  ST  AT  I C  PRESSURE  AT  1ST  STAGE  INLET  (BARS) 

VC  10) =TOT AL  PRESSURE  AT  1ST  STAGE  INLET ( BARS  ) 

(BRITISH  UNIT S*PSI A ) 

(8RITISH  UNIT S=P5I A  ) 

VCI1)=STATIC  TEMPERATURE  AT  1ST  STAGE  INLETCDEG  K) 

(BRITISH  UNIT S=OEG  R) 

VC  12 )  =  TliT AL  TEMPERATURE  AT  1ST  STAGE  INLET  (DEG  K) 

(BRITISH  UNIT S=OFG  RI 

V()3)*GAS  DENSITY  AT  1ST  STAGE  INLET  CKG/CU.M) 

(BRITISH  UNITS*LBM/CU.FT  ) 

V(14)rGAS  VELOCITY  AT  1ST  STAGE  NOZZLE  EXIT  CM/SEC) 

(BRITISH  UNITS  =  FT/SEC  * 

V C 1 5 )  =  MEAN  SPECIFIC  HEAT  OF  GAS  THROUGH  1ST  ST AGEC JOULES/KG-DEG  K) 
(BRITISH  UNITS=BTU/LRM-OEG  R) 

V ( 16 ) =MF AN  SPECIFIC  HEAT  RATIO  THROUGH  1ST  STAGE 
VC 17)«S0DERBURG  LOSS  COEFFICIENT  FOR  1ST  STAGE  ROTOR 
VC I8)=S0DERBURG  LOSS  COEFFICIENT  FOR  1ST  STAGE  STATOR 
VC  19) =TOT AL  TO  TOTAL  EFFICIENCY  FOR  1ST  STAGE 
VC  20) x AX  I AL  VELOCITY  OF  GAS  IN  1ST  STAGE  (M/SEC) 

(BRITISH  UNITS=FT/SEC  ) 

VC  21 ) *  ST  AT  j  C  PRESSURE  AT  1ST  STAGE  ROTOR  EXIT  (BARS) 

(BRITISH  UNITS  =  PSI A  ) 

VC  22 ) sTOTA  L  PRESSURE  AFTER  1ST  ST AGE ( BAR S  )  ( BR I T I SH  UNITS-PSIA) 
V(23)»STATIC  TEMP.  AT  1ST  STAGE  ROTOR  EXIT  (DEG  K) 

(BRITISH  UNITS=DFC,  K) 

VC  24) =T  OT  AL  TEMP.  AFTER  1ST  STAGE  (DEG  K )  (BRITISH  UNITS*DEG  R) 

VC  25 ) =GAS  DENSITY  AT  1ST  STAGE  ROTOR  EXIT  KG/CU.M) 

(BRITISH  UNITS  =  Lf'M/CU.FT  ) 

VC  26 )  =  GAS  VFLOCITY  AT  1ST  STAGE  ROTOR  EXIT  (M/SEC) 

(BRITISH  UNITS=FT/SFC ) 

V C 27 )  =  ANGL E  OR  RELATIVE  f,AS  VELOCITY  VECTOR  AT  1ST  STAGE  NOZZLE 
EXIT  (RADIANS)  (BRITISH  UN  I TS“DEGRE ES ) 

VC  28  )  =  ANGL  F  OF  ACTUAL  GAS  VELOCITY  VECTOR  AT  1ST  STAGE  ROTOR  EXIT 
CRAOIANS)  (BRITISH  'JNI TS-OEGREES ) 

VC 29)=EFFFCTIVE  ROTOR  EXIT  AREA  OF  1ST  STAGE  CSO.M) 

(BRITISH  UNIT S=SO.FT ) 

V(30)*P0WER  OUTPUT  OF  1ST  STAGE  (MW)CBRITISH  UNI TS*HORSEPOWER ) 
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VI 31) “STATIC  PRESSURE  AT  2ND  STAGE  1NLFT  (BARS) 


SECOND  STAGE 

IBRITISM  UNIT S-PSI A  ) 

VI32)«ST4T1C  TEMPERATURE  iT  2ND  STAGE  INLETIOEG  Kl 
•BRITISH  UNITS-DEG  RI 

VI  331-GAS  DENSITY  AT  2ND  STAGE  INLET  IKG/CU.M) 

•  BRITISH  UNITS-LBM/C.  .FT  ) 

VI  341-GAS  VELOCITY  AT  2ND  STAGE  NOZZLE  EXIT  IH/SEC) 

•British  un:ts=et/sec) 

VI  351-MEAN  SPECIFIC  HEAT  OF  GAS  THROUGH  2ND  STAGE  I JOULES/KG-DEG  K> 
•BRITISH  UNITS-BTU/LRM-DEG  R) 

V 1 341-MEAN  SPECIFIC  HEAT  RATIO  THROUGH  2ND  STAGE 
VI37)«SO0FRRURG  LOSS  COEFFICIENT  FOR  2ND  STAGE  ROTOR 
V 1 38 )- SOOERRtlRG  LOSS  COEFFICIENT  FOR  2ND  STAGE  STATOR 
VI 39) -TOTAL  TO  TOTAL  EFFICIENCY  FOR  2ND  STAGE 
VI 40 ) -AXIAL  VELOCITY  OF  GAS  IN  2ND  STAGE  IM/SEC) 

•BRITISH  UNITS-FT/SEC) 

VI41)- STATIC  PRESSURE  AT  2N0  STAGE  ROTOR  EXIT  IBARS) 

•  BRITISH  UNIT S-PSI A  ) 

VI42)«T0TAL  PRESSURE  AFTER  2ND  STAGE  I  BARS) IBRI TI SH  UNITS-PSIA) 

VI 43) -STATI C  TEMP.  AT  2ND  STAGE  ROTOR  EXIT  IDEG  K> 

•  BRITISH  UNITS-DEG  R  ) 

VI 44)*T0TAL  TEMP.  AFTER  2ND  STAGE  IDEG  K)  IBRITISH  UNITS-DEG  *! 
VI45)«GAS  DENSITY  AT  2ND  STAGE  ROTOR  EXIT  KG/CU.M) 

IBRITISH  UNI TS-LBM/CU.FT  ) 

VI 46 )*GAS  VELOCITY  AT  2ND  STAGE  ROTOR  EXIT  IM/SEC) 

IBRITISH  UNITS=FT/SEC> 

VI47 )- ANGLE  OF  RELATIVE  GAS  VELOCITY  VECTOR  AT  2ND  .STAGE  NOZZLE 
EXIT  IRADIAMS)  IBRITISH  UNI TS-DEGREES ) 

V 1 48 )- ANGLE  OF  ACTUAL  GAS  VELOCITY  VECTOR  AT  2ND  STAGE  ROTOR  EXIT 
•RADIANS)  IBRITISH  UNI TS-DEGREES > 

VI 49 INEFFECTIVE  ROTOR  EXIT  AREA  OF  2ND  STAGE  ISO. HI 
IBRITISH  UNIT S-SO.FT  ) 

V 1 501-POWER  OUTPUT  OF  2ND  STAGE  IMWMBRITISH  UNI  TS-HORSEPOWER  ) 
THIRD  STAGE 

VI5II-STATIC  PRESSURE  AT  3RD  STAGE  INLET  IBARS) 

IBRITISH  UNITS-PSIA) 

VI52I-STATIC  TEMPERATURE  AT  3RD  STAGE  INtETIDEG  K) 

IBRITISH  UNITS-DEG  R) 

VI53I-GAS  DENSITY  AT  3RD  STAGE  INLET  (KG/CU.M) 

IBRITISH  UNITS-LBM/CU.FT I 

VI54J-GAS  VELOCITY  AT  3RD  STAGE  NOZZLE  EXIT  (M/SEC) 

IBRITISH  UNITS-FT/SEC) 

VI 55 ) -MEAN  SPECIFIC  HEAT  OF  GAS  THROUGH  3RD  STAGEI JOULES/KG-OEG  K> 
IBRITISH  UNITS-BTU/LSM-OEG  R> 

VI56)-MEAN  SPECIFIC  HEAT  RATIO  THROUGH  3RD  STAGE 
VI 57) -SOOFRbURG  LOSS  COEFFICIENT  FOR  3RD  STAGE  ROTOR 
VI 58 )=SODERBURG  LOSS  COEFFICIENT  FOR  3RD  STAGE  STATOR 
VI 59)  =  T OTAL  TO  TOTAL  EFFICIENCY  FOR  3RD  STAGE 
VI 60) -AXIAL  VELOCITY  OF  GAS  IN  3RD  STAGE  (M/SEC) 

•BRITISH  UNITS-FT/SEC) 

VI 61 l-STAT IC  PRESSURE  AT  3RD  STAGE  ROTOR  EXIT  (BARS) 

VI 62) -TOTAL  PRESSURE  AFTER  3RD  STAGE  I  BAR S ) < BRJ TI SH  UNITS-PSIA) 
IBRITISH  UNITS-PSIA) 

VI 63)*STATIC  TEMP.  AT  3RD  STAGE  ROTOR  EXIT  IDEG  K) 

IBRITISH  UNITS-DEG  R) 
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V(64)«T0TAL  TEMP.  AFTER  3P0  STAGE  I  DEG  K )  IBRITISH  UNITS«DEG  R) 
V(65)*GAS  DENSITY  AT  3RD  STAGE  ROTOR  EXIT  KG/CU.M) 

(BRITISH  UNITSsLEM/CU.FT  ) 

V(66)»GAS  VELOCITY  AT  3RD  STAGE  ROTOR  EXIT  (M/SEC) 

(BRITISH  UN!TS=FT/SEC) 

V 1 67 1 * ANGl F  OF  RELATIVE  GAS  VELOCITY  VECTOR  AT  3RD  STAGE  NOZZLE 
EXIT  (RADIANS  I  (BRITISH  UN  I TS  =  DFGREES  ) 

V(  68 1  “ANGL E  OF  ACTUAL  GAS  VELOCITY  VECTOR  AT  3RD  STAGE  ROTOR  EXIT 
(RADIANS)  (BRITISH  UNI «S“DEGREES ) 

V( 69 (“EFFECTIVE  ROTOR  EXIT  AREA  OF  3RO  STAGE  (SO.H) 

(BRITISH  UNI TS“SD.FT ) 

V(70)«P0WFR  OUTPUT  OF  3RD  STAGE  (MWXBRITISH  UNI  TS=HORSEPOWER ) 


EXHAUST  SYSTEM  VARIABLES 


V(71)» STATIC  PRESSURE  IN  EXHAUST  DUCT  IBARS) 

(BRITISH  UNI TS=PSI A  ) 

V(72)* STATIC  TEMP.  IN  EXHAUST  DUCT  (DEG  KMBRITISH  UNI  TS“DEG  R) 

V( 73) “GAS  DENSITY  IN  EXHAUST  DUCT  (KG/CU.M) 

(BRITISH  UNITS“LBM/CU.FT) 

V( 74) =GAS  VELOCITY  IN  EXHAUST  OUCT  (M/SEC) 

(BRITISH  UNITS=FT /SEC) 

V(75)-STATIC  PRESSURE  AT  EXHAUST  NOZZLE  THROAT  (BARS) 

(BRITISH  UNITS=PSI A) 

V( 76) “STATI C  TEMP.  AT  EXHAUST  NOZZLE  THROATt DEG  K) 

(BRITISH  UNITS=DEG  R) 

V( 77) =GAS  DENSITY  AT  EXHAUST  NOZZLE  THROAT  (KG/CU.M) 

(BRITISH  UNITS=LBM/CU.FT) 

V( 78) “GAS  VELOCITY  AT  EXHAUST  NOZZLE  THRDAT  (M/SEC) 

(BRITISH  UNITS=FT /SEC) 

V(79)=MEAN  SPECIFIC  HFAT  OF  EXHAUST  GAS  (JOULES/KG-DEG  K) 
(BRITISH  UNITS  =  BT'J/LBM-DEG  R) 

V ( BO) “MFAN  SPECIFIC  HEAT  RATIO  OF  EXHAUST  GAS 

V(Bl)«EFFECTIVE  EXHAUST  NOZZLE  THROAT  AREA  (SO.M) 

(BRITISH  UNITS=SQ.FT ) 
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c 


0001 

DIMENSION  VI8 1 t*R (81 > .PDIOl. 81 I.VCORR 181 I.OESI 81 1.0(811 

0002 

COMMON  A.Z.O .ALF2.ALF4 .ALF&.ALFR , ALF10  « BET3 «BET5« BET7.8ET9.6ETH 

0003 

COMMON  HOB  1 .H3fi2.H0B3.H0B4, H9B5 .Ul.U2.U3.U4.05 »A 1 . A3 «A5 <A7 , A9 

0004 

COMMON  AT .ADUCT , ALF1 

000$ 

c 

COMMON  ALTO.GftBE • ALTE 

e 

c 

READING  THE  OATA  CAROS 

0006 

RE AO ( 5 , 10 )  NV Aft 

0007 

10 

FORMAT (131 

0008 

READ (5.11)  TLRNCE.  ITMAX 

0009 

11 

FORMAT (F 10.4 ,  110) 

0010 

RE AD  1 5 . 12 )  (V  ( I  )  •  I  =  l.NVAR) 

0011 

12 

FORMAT (BF 10.3) 

0012 

REAO  (5.13)  (DES(I).  I  =»  1.  NVAR ) 

0013 

13 

FORMAT  ( 20A4  ) 

0014 

00  1000  LM*1»  5 

001$ 

ALTO=2.0-.4*(LM-1) 

0016 

ALTE* .9 

0017 

GRBE* .95 

0018 

D*.254 

0019 

2*. 203 

0020 

A*. 003048 

0021 

ALF1*0. 

0022 

ALF2*1.29 

0023 

BET3*=1  .03 

0024 

Ul-516. 

002$ 

Al*. 000229 

0026 

H0B1=1.25 

0027 

ALF4* 1.2277 

0026 

BET5M.2277 

0029 

U2*S16. 

0030 

A 3*. 000567 

0031 

H0B2=l .25 

0032 

ALF6-1.2277 

0033 

B6T7=1.2277 

0034 

AS*. 00102 

003$ 

U3-516. 

0036 

H0B3-1.25 

0037 

ADUCT=. 00811 

L 

c 

r 

WRITING  OUT  THE  INPUT  OATA 

0038 

WRITE (6,20)  NVAR 

0039 

20 

FORMAT ( 1H1 ,///»•  NUMBER  OF  VARIABLES  *  14) 

0040 

WRITE(6,21)  TLRNCE 

0041 

21 

FORMAT ( '  MAXIMUM  FRACTION  CHANGE  FOR  CONVERGENCE  =  '.  F10.4,  It) 

0042 

WRI TE (6 , 22 ) 

0043 

22 

FORMAT ( 'OVARI  ABLE  NUMBER  AND  ITS  TRIAL  VALUE') 

0044 

WRITE  (6,23)  (J,  DESU).  V(J),  J  =  1,  NVAR) 

004$ 

c 

23 

FORMAT  ('  V ( * .  13.  ')  *  A4»  •  *  F15.5) 
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0046 


0047 


0048 


0049 

0050 

0051 


0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 

0061 

0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 

0081 


C  INITIALIZING  THE  ITERATION  COUNTER 
ITER  *  1 

C  CALLING  SUBROUTINES  TO  CALCULATE  VALUES  OF  RESIDUALS*  PARTIAL 
C  DERIVATIVES  AND  CHANGES  IN  VALUES  OF  VARIABLES 

C 

30  CALL  EONS (NVAR,  V.  RJ 
C  WRITE  (6.33  ) 

C  33  FORMAT! 'OEOUATION  NUMBER  RESIOUAL* ) 

C  WRITE  (6,35)  (I  ,  R(I).  I  «  l.NVAR) 

C  35  FORMAT  (110.  F20.5) 

CALL  PAPQ IF (NVAR ,  V,  ft,  PO) 

C 

C  PRINTING  OUT  NON-ZERO  VALUES  OF  PARTIAL  DERIVATIVES 
C 

DO  54  1=1.  NVAR 

DO  53  J  =  1,  NVAR 
ZT=  ABSlPOd  .J)  ) 

C  REMOVE  THE  C  FROM  THE  NEXT  THREE  CAROS  IF  PRINTOUT  DESIREO 

C  IF ( ZT-  O.COOOOOOl)  53,  53,  51 

C  51  WR I TE (6 , 52 )  I.  J.  PD(I.J) 

C  52  FORMAT ( •  PD(',I2,','.I2.')  =•  ,  F20.10) 

53  CONTINUE 

54  CONTINUE 

CALL  GAUSSY(PO.  R,  VCORR.  NVAR) 

C  CORRECTING  THE  VALUES  OF  THE  VARIABLES 
DO  44  L  =  l.NVAR 
44  VIL)  =  V ( L )  -  VCORR(L) 

C 

0(  1)«V(  1J/745.7E-06 
Q(  2  I  =V (  2  1*2.20462 
0(  3  )=V(  3)*2. 20462*. 092903 
Q(  4  )=V  (  4) 

0<  5 )=V (  5) 

0(  6 )=V i  6) 

0(  7 )=V (  7) 

0(  8 )=V (  8  ) 

0(  9  )=V  (  91/6894. 757E-05 
0( 10)=V( 10J/6694.757E-05 
0( 1 1  )=V ( 111*1.8 
0(12  )=V ( 12  )*1 .8 
0( 1 3  )=V ( 131/16.02 
0(14)  =V (14)/. 3048 
0(15  )=V(  155  /4184. 

0(16  )=V( 16 ) 

0(17  )=V (17) 

0(1©)  = V (18) 

0(19  )=V ( 19) 

0(20)=V(20)/.3C48 
0(21  )=V( 211 /6894.757E-05 
0(22)=V(22 1/6894 . 757E-05 
0( 23)=V( 23)*1.8 
0(24  )=V(24  )*1 .8 
0(25)=V(25)/16.02 
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0062 

0063 

0064 

0065 

0066 

0067 

0086 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0096 

0099 

0100 

0101 

0102 

0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 

0111 

0112 

0113 

0114 

0115 

0116 

0117 

0118 

0119 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 


0«26)=V(26)/.3048 
0(27l«V(27  1*57.3 
0(28)=V(28**37.3 
0(29)=v(291/«. 3048**2) 
0(30»«V (301/745. 76-06 
0(3l)sV«31 1/6894.7576-05 
0( 32  )*V( 32 1  *  1, 8 
0(33  )*V( 33 1/16.02 
0(341=V(34J/.304P, 

0 (35 1*V « 35  1/4184. 
0(361=V«36J 
0(37 1*V ( 37  1 
0(381*V(38) 

0(39 J=V«39) 

0(401=V«40)/.3048 
0(411=v«41  1/6894.7576-05 
0(42)*V(42 1/6394. 7576-05 
0(43  1*V(43 1*1.8 
0 ( 44 1 * V ( 44 ) *  1 , 8 
0(451=V(45 1/16.02 

0(46I=V(46) /.304? 

0(47  )*V(47  1*57.3 
0(481=V(48)*57.3 
0(49|=v<’49)/(  .3048**2) 

0  501=V«50)/745.7E-06 
0 ( 52  1=V ( 52  1*1.8 

°tf!J=vt511/6894.757E-05 

0(53  1»V (53 1/16.02 
0(541*V(541/.3048 
0 ( 55  1*V (55  1/4184. 

0( 56  )=V( 56 ) 

0 (57 1 =V (57  1 
0( 58 )=V ( 5  3  1 
0 (59  1  *V ( 59 ) 
0(60J=V(60)/.3048 
0(61)  =v (61  1/6  894.7576-05 
0 ( 62 1=V( 62 1 /6894. 7576-05 
0 (63 1=V (63  1*1,8 
0 (64 |sv { 64  ) *  1. 8 
0(65 1*V (65 1/16.02 
0(661=V(66)/,3048 
0 (67  )  =V (67  1*57.3 
0(68J=V(6fi)*57.3 
0(69 )=V(69 )/« .3048**2) 
0(70)*V«70)/745. 76-06 
(71  I  =V (71  1/6394.7576-05 
0«  72 1  =  V( 721*1.8 
0(73  )=V(73 1/16.02 
01 74 )=V ( 74 1/.3048 
0(75  )  =V ( 75 J/6394.75  7f_05 

0«76)=V(76)*1.8 

0(77  l=v (77  1/16.02 
0(781=V(78)/.3048 
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0135  0(79)*V(79)/4184. 

0136  0(80)=V(80> 

0137  0(81)-V(81)/(. 3048**2) 

C 

C  WRITING  OUT  RESULTS  OF  THIS  ITERATION 
C 

C  TO  PRINT  OUT  EVERY  ITERATION  REMOVE  COMMENT  NOTATION  FROM 
C  FOLLOWING  CARDS 

C 

C  WRITING  OUT  RESULTS  OF  THIS  ITERATION 
C 

1  C  WRITE  16.321  ITER 

C  32  FORMAT  PORESULTS  AFTER*.  14.  *  ITERATIONS') 

C  WRITE (6*34) 

C  34  FORMATI • O', 4X,» VAR! ABLE' .14X,  'VALUE ».9X. 'CHANGE  FROM  PREVIOUS' » 

C  WRITE  (6.36)  (I.  DESU),  V(I>.  VCORR(I),  I  *  l.NVAR) 

C  36  FORMAT  ('  VI'.  12.  ')  =  A4,'  «  2F21.5) 

C 

C  TERMINATING  IF  MAXIMUM  NUMBER  OF  ITERATIONS  REACHED  OR  OTHERWISE 
C  INCREMENTING  THE  ITERATION  COUNTER 

0138  IF ( ITER  -  IT MAX )  3R,  99,  99 

0139  38  ITER  =  ITER  *  1 

C 

C  CHECK  TO  SEE  IF  CHANGE  OF  VARIABLE  IS  LESS  THAN  SPECIFIED  TOLERANCE 
C 

0140  K  •  1 


0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 


40  VAL  *  ABSIVCORRIK))  -  ABS! TLRNCE*V( K) ) 

IF(VAL)  41,  30.  30 

41  IF  IK  -  NVAR)  42,  99,  99 

42  K  *  K  ♦  1 
GO  TO  40 

99  WRITE (6, 122) 

1 22  FORMATI 1H1) 

WRITE (6,111) 

111  FORMAT  COVAR I  ABLE  NUMBER  AND  ITS  FINAL  VALUE'/) 

WRITE  16, 1973 )  ALTO 

1973  FORMAT! •  THESE  VALUES  ARE  FOR  AN  ALTERNATOR  OUTPUT  OF  '.F4.2,'  MW' 
1/) 

WRITE (6, 1 12  > 

112  FORMATI 'O' ,4X , *VARI ABLE ', 12X, ' VALUE! SI ) • ,21X,' VALUE! BRIT ISH ) • ) 

WR ITEI6, 113)1 NM,OES!NH),V(NM), 01 NM) , NM= 1 , 15  ) 


FORMAT!'  VC 
IP'./. 

.12.  •)  = 

• 

,  A4,«  «  • 

.F14.5, 

•  MW 

•  .F14.5, '  H 

2* 

V! 

'.12,') 

=  ' , A4 , ' 

S 

• ,F 14.5.  ’ 

KG/ SEC 

•  .6X.F14.5,' 

LBM/SEC  './. 

3* 

40 

V! 

.FT 

• .1 2.  • ) 

' » / » 

=  •  ,  A4  i  ' 

K 

• ,F 14.5, ' 

KG/SEC- 

SO. M* .2X.F14 

.5.'  LBM/SEC-S 

5' 

V! 

'.12,') 

=  • ,  A4  »  ' 

« 

• »F14. 5, • 

* ,6X, F14.5, ' 

'./♦ 

6* 

VC 

'.12.'  ) 

=  ' , A4 , ' 

s 

• .F14.5. • 

».6X,F14.5.» 

•  ♦/. 

7* 

V! 

'.12.  • ) 

=  ' ,  A4  » ' 

X 

• .F14.5. • 

'.6X.F14.5, • 

'./, 

8* 

V! 

'.12.'  ) 

=  ' , A4, ' 

s 

• ,F14.5» ' 

•  .6X.F14.5,' 

',/. 

9' 

V! 

'.12.  •) 

s  '  ,  A4 ,  ' 

s 

•  ,F14.5, ' 

KG/MOLE 

• ,6X, F14.5, ' 

LBM/MOLE*.//, 

1' 

V! 

• .12,'  ) 

«  * , A4  » ' 

s 

• .F14.5. ' 

bars 

• .6X.F14.5,' 

PSIA  ',/, 

2* 

V! 

'.12,') 

=  '  ,  A4  ,  ' 

X 

• ,F 14. 5,  • 

BARS 

' ,6X, F14.5, ' 

PSIA  ',/, 
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3*  Vi*.I2.*>  *  «.A4,'  = 
4*  VI  '  ,12, ' )  *  • ♦ A4« '  = 
5*  VI ' *12* '  I  *  • « A4 ,  •  « 
6*  VI*.  12.'  )  *  • «  A4 , •  = 
7*  VI*. 12.*)  =  • • A4 • •  = 
8G  RM 

WRITE  16.1 14 MNM.OESINM) 


FORMAT! • 

VI 

*  9 

12.  *) 

*  *,  A4,»  = 

•  .F14.5, 

•  *F  14 

.5 

./. 

1' 

VI 

•  .12* 

•  ) 

S 

* .  A4 ,  • 

=  ' ,F 14. 5  .  * 

'.6X.F14.5.  • 

./, 

2* 

VI 

•  .12. 

•  ) 

= 

'  . A4 ,  ' 

=  *,F14.5.* 

•  .  6X  .  F14 . 5  *  • 

•  /  « 

3* 

VI 

*.I2. 

• ) 

s 

•  .  A4  .  * 

=  * .F 14. 5  .  ' 

'.6X.F14.5.' 

./, 

4* 

VI 

*  *  1 2, 

= 

’  . A4 ,  • 

=  *  »F 14,5. * 

M/SEC 

• .6X.F14.5.' 

FT/SEC 

./ 

5* 

VI 

•  .12. 

•  i 

= 

• . A4 .  • 

=  '.F14.5.' 

BARS 

•  .6 X.  FI  4. 5.  • 

PSIA 

•  / » 

6* 

VI 

•  .12. 

• » 

= 

'.A4.' 

=  *.F14.5.' 

bARS 

• .6X.F14.5.' 

PSIA 

»/ » 

7* 

VI 

•  .12  ♦ 

•  > 

= 

•  «  A4 .  • 

*  • .F 14.5. * 

DEG  K 

• *6X*  F14.5. * 

DEG  R 

./ 

8* 

VI 

•  .  1 2. 

• ) 

s 

•  . A4 ,  • 

=  •  .F 14.5. * 

DEG  K 

• .6X.F14.5.' 

DEG  R 

./ 

9* 

VI 

•♦12. 

' ) 

= 

•  «  A4 ,  • 

*  * .F 14. 5 . * 

KG/CU.M 

• .6X. F14.5. * 

LBM/CU. 

FT 

•  ./ 

1* 

VI 

•  .12. 

•  > 

s 

'.  A4,* 

*  '.F14.5.* 

M/SEC 

• .6X.F14.5.' 

FT/SEC 

./ 

2* 

VI 

•  .1 2. 

•  > 

T 

•  »  A4  «  • 

=  *  »F 14. 5. ' 

RADIANS 

* »6X. F14. 5. * 

DEGREES 

./, 

3' 

VI 

•♦12. 

•  > 

S 

•  ,A4,' 

=  *.F14.5.' 

RADIANS 

•♦6X.F14.5.' 

DEGREES 

./. 

4* 

VI 

•  .12. 

•  > 

S 

•  .  A4  ,  • 

*  *  *F 14. 5 « * 

SO  .M 

'.6X.F14.5. * 

SO. FT 

./, 

5* 

VI 

'.12. 

•  > 

S 

•  .A4,* 

=  *.rl4.5.' 

MW 

• .6X.F14.5.' 

HP  • 

/) 

WRITE (6. 115 MNM.OESINM) 


*fl4.5.'  DEG  K  • *6X*  F14.5*  *  DEG  R  *./, 
*F 14.5. •  DEG  K  *.6X,Fi4.5.'  DEG  R  './. 

* F 14. 5* •  KG/CU.M* .6X.F14.5,'  LBM/CU .FT * ,/ * 
* F 14. 5. •  M/SEC  • , 6X ,F 14 .5 • '  FT/SEC  •  ,/. 
* F 14. 5 * •  J/KG-DEG  K*, 3X.F14.5*'  BTU/LBM-DE 

VINM)  .OINM)  ,NM*16.30) 


V I NM ) .0 1 NM ) , NM=  3 1 


115  FORMAT! •  VI 

*  f 

12.  •) 

5  • 

.  A4,  •  * 

*  .F14.5 

.'  BARS 

•  ,F14 

.5 

.•  p 

1SIA* 

./. 

2* 

VI 

',12,'  ) 

S 

' , A4 , ' 

s 

* ,F14.5,  • 

DEG  K 

'  ,  6X  ,  F14 . 5  ,  • 

DEG  R 

• 

» /  • 

3* 

V! 

• .1 2.  •  1 

= 

• , A4 ,  • 

= 

* ,F 14. 5 , ' 

KG/CU. 

M',6X,F14.5,' 

LBM/CU. 

FT 

• ,/. 

4* 

VI 

*  .12.  •  I 

= 

' , A4 , ' 

s 

' ,F 14.5,  • 

M/SEC 

•  .6X.F14.5,' 

FT/SEC 

• 

./. 

5* 

VI 

• ,1 2 , ' ) 

= 

' , A4 ,  • 

s 

• ,F 14. 5, • 

J/KG-DEG  K'.3X,F14. 

5,'  BTU/LBM-DE 

6G 

R* 

./. 

7* 

VI 

* .  1 2 .  • ) 

s 

•  .  A4 ,  • 

9 

*,F14.5,» 

• ,6X, F14. 5, ' 

.  /» 

8* 

VI 

'.12.* ) 

s 

• ,  A4, ' 

= 

' ,F 14.5, ' 

• .6X.F14.5,' 

./. 

9* 

VI 

*.I2,*» 

= 

•  , A4 ,  • 

S 

•  ,F 14. 5 , ' 

'.6X.F14.5,' 

V  /» 

1* 

VI 

'.12,'  ) 

s 

•  ,  A4, ' 

E 

»,F14.5,  • 

'.6X.F14.5,' 

,/. 

2* 

VI 

'.12,  ') 

s 

• , A4,  ' 

= 

' ,Fl4.5, ' 

M/SEC 

'  ,6 X, F14.5,  ' 

FT/SEC 

) 

WRITEI6, 11611 HM.OESINM) 

,VINK),OINM),NM= 

41,50  ) 

116  FORMAT  I*  VI 

•  t 

12.  •) 

*  • 

,  A4 »  '  a 

'  .F14.5 

,'  BARS 

•  ,F14 

.5 

.  •  ■' 

1SIA* 

./. 

2* 

VI 

'.12,  •) 

; 

• , A4 ,  • 

S 

» ,F  14. 5  • » 

BARS 

'  ,6X, F14.5,  ' 

PSIA 

• 

./. 

4* 

VI 

'.12.'  ) 

= 

' , A4,  • 

= 

' ,F14.5.  ' 

DEG  K 

•.6X.F14.5,' 

DEG  R 

• 

»  /  , 

5' 

VI 

'.12.  'I 

= 

•  , A4  ,  • 

= 

' »F 14. 5.  • 

DEG  K 

'.6X.F14.5,  • 

DEG  R 

• 

,/, 

6* 

VI 

',12,'  » 

= 

', A4.' 

= 

'.F14.5. ' 

KG/CU. 

M* .6X.F14.5 ,' 

LBM/CU. 

FT 

•./, 

7* 

VI 

'.12,  ') 

= 

' , A4 , ' 

S 

« ,F 14. 5, ’ 

M/SEC 

'.6X.F14.5.' 

FT/SEC 

• 

./, 

8* 

VI 

',12.'  ) 

X 

' ,  A4 , ' 

= 

» ,F 14, 5, ' 

RADI  AMS' , 6X.F14.5,' 

DEGREES 

• 

,/. 

9* 

VI 

' .  12  , '  1 

= 

•  , A4  ,  • 

= 

'.F14.5.' 

RADI ANS* ,6X, F14.5. • 

DEGREES 

• 

,/. 

1* 

VI 

'.12.'  » 

E 

•  . A4,' 

E 

'.F14.5.' 

SO.M 

' ,  6X  ,F  14 . 5 . ' 

SQ.FT 

• 

»  /  » 

2* 

VI 

',12,') 

= 

'  , A4 , ' 

= 

• ,F 14.5, • 

MW 

' »6X, F14.5, • 

HP  ' 

9 

/> 

WRITE16.117)INM,DES(NM| 

, V  I NM ) , 0 ( NM ) , NM= 

51,60  ) 

117  FORMAT! •  VI 

1  9 

12,  '  ) 

s  1 

,  A4.'  = 

'.F14.5 

,  •  B  AR  S 

•  ,F14 

.5 

.•  p 

1SIA ' 

./. 

2' 

VI 

',12.  •) 

= 

* , A4 ,  • 

= 

* ,F 14. 5, ' 

DEG  K 

'  ,6 X,  FI 4.  5, ' 

DEG  R 

• 

,/. 

3* 

VI 

'.12.') 

9 

•  ,A4.  • 

- 

'.F14.5. ' 

KG/CU. 

M'  ,  6X , F 14 . 5  . ' 

LBM/CU. 

FT 

*./. 

4* 

VI 

',12,') 

S 

• . A4 ,  • 

s 

’ ,F 14. 5,  ' 

M/SEC 

•.6X.F14.5,' 

FT/SEC 

• 

./. 
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5*  Vl»»12,'> 
6G  R*»/i 

X 

• .  A4,  • 

9 

' ,F 14. 5 ,  • 

J/KG-DEG  K • , 3X  ,F14. 

5,'  BTU/LBM-DE 

7*  VI*. 12.*) 

X 

• , A4 ,  • 

9 

• ,F 14. 5, • 

• »6X» FI  4.5. • 

8*  Vt'»!2»'> 

9 

• , A4. • 

9 

'.F14.5,' 

• ,6X,F14,5,' 

9*  VI'. 12.') 

9 

•  . A4 , ' 

9 

'.F14.5.' 

• ,6X. FI  4, 5 . ' 

l*  VI '.12, ' > 

9 

•,A4,' 

9 

• .F14.5. • 

• .6X.F14.5.' 

2*  VI • « 12* • l 

9 

• ,  A4 ,  • 

r 

'.F14.5,' 

M/SEC  1 .6 X,  FI  4. 5,  • 

FT/SEC  '  ) 

0164 

WRITE  16. 118  UNM.DESINMl 

, V 1 NM ) . 0 ( NM ) , NM=61 . 70 ) 

0165 

118  FORMAT ( •  V( 
1S1A' »/» 

•t 

12,  •> 

X  • 

,  A4.  •  * 

•  ,F 14. 5. '  BARS 

• »F14.5«  '  P 

2*  V('.I2.') 

= 

•  , A4  ,  • 

s 

•  ,F 14. 5 ,  • 

BARS  ' .6 X, F 1 4, 5, ' 

PSIA  './, 

4*  V( '.I?.'  ) 

S 

*  .  A4  ,  • 

= 

'.814.5,  • 

DEG  K  '.6X.F14.5.' 

DEG  R  •  ./. 

5*  VI'.I?,') 

X 

•  ,A4,  • 

X 

• »F 14. 5.  • 

DEG  K  '.6X.F14.5,' 

DEG  R  '  ./. 

6*  V I • . I 2*  •  » 

X 

•,A4,' 

X 

' ,F 14.5,  • 

KG/CU.M'  .6X.F14. 5,' 

LBM/CU.FT • ,  / . 

7*  VI • .12,  •  ) 

X 

•  .  A4  ,  ' 

s 

' ,F 14. 5  .  • 

M/SEC  '.6X.F14.5.' 

FT/SFC  •  ,/. 

8*  VI  *  » 1 2.  •  ) 

X 

'  , A4,' 

r 

' ,F14«5, • 

RADI ANS* .6X.F14.5,' 

DEGREES  './. 

9*  VI ' ,12,  ‘ ) 

X 

•  ,  A4 , « 

x 

' »F14.5.  • 

RADIANS', 6X.F14.5,' 

DEGREES  './, 

1*  VI '.12,'  ) 

= 

•  ,A4,' 

X 

• ,F14.5.  • 

SC'.M  '»6X,F14.5»' 

SO. FT  ',/. 

2*  VI*. 12,') 

= 

•  , A4  ,  • 

X 

• ,F 14. 5 ,  • 

MW  ' «6X. FI  4. 5, ' 

HP  •  ,/) 

0166 

WRITEI6»119)INM»DESINM) 

,  V I NM ) , 0 1 NM ) , NM=  71,81  ) 

0167 

119  FORMAT!'  VI 
1SIA'./. 

12.  •) 

X  • 

,  A4, '  = 

• ,  F 14. 5, •  BARS 

' ,F14.5, •  P 

2'  VI*. 12.') 

X 

• . A4  ,  ' 

x 

' ,F 14. 5.  • 

DEG  K  • ,6 X, F 1 4.5, ' 

DEG  R  •  ,/, 

3*  VI ' , 1 2, '  ) 

s 

• . A4 ,  • 

X 

•  »F14.5. • 

KG/ CU.M* , 6X.F14.5,' 

LBM/CU.FT',/, 

4*  VI' ,12. ') 

X 

'  ,  A4  ,  • 

= 

• ,F 14, 5,  ' 

M/ SEC  '  ,6X, FI  4. 5, ' 

FT/SEC  '  ,/, 

5*  VI', 12,'  ) 

X 

',  A4,  ' 

= 

• ,F 14.5.  ' 

BARS  • ,6X ,F 14 .5 , ' 

PSIA  ',/. 

6'  VI*. 12,') 

= 

«  .  A4  ,  • 

X 

• ,F1 4. 5 . ' 

DEG  K  '.6X.F14.5, ' 

DEG  R  •  ,/. 

7*  VI  ',12, '  ) 

x 

• , A4  ,  ' 

= 

• .F14.5,  • 

KG/CU.M'. 6\,F14. 5,' 

LBM/CU.FT ',/. 

8*  VI', 12,') 

x 

' , A4  , ' 

X 

• ,F 14. 5. » 

M/SEC  ' ,6  X. FI  4.5, * 

FT/SEC  '  ,/, 

9*  VI', 12.' ) 
1G  R'./. 

'  »  A4  ,  ' 

s 

• , F 14, 5,  ' 

J/KG-DEG  K' ,3X, F14. 

5,'  BTU/LBM-DE 

2*  VI', 13,'  ) 

= 

•  ,  A4  ,  ' 

= 

' »F 14. 5,  • 

'.,  6X,F  14 . 5 , ' 

'./, 

0168 

0169 

0170 

3'  VI ',13.  •) 
1000  CONTINUE 
STOP 

END 

•  , A4  ,  ' 

• ,F 14. 5, • 

SO .M  ',6X,ri4,5,' 

SO, FT  •) 
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0001 

0002 


0003 


0004 

0005 

0006 

0007 

0000 

0009 

0010 

0011 

0012 

0013 

0014 

0015 


SUBROUTINE  PARDJF  ( NVAR ,  V,  R,  PO) 

DIMENSION  V(Bl),Riei).PD<ei.81),VDIBl).RD(Bl» 

c 

C  GLOSSARY  FOR  SUBROUTINE  PAROIF 

C  OV  =  FRACTION  OF  VARIABLE  CHANGE  USED  IN  TAKING  PARTIAL 

C  vr>  «  V  +  DELTA  V  «  V  ♦  V*DV 

C  RD  *  R  EVALUATED  AT  VD 

C  THE  PARTIAL  DERIVATIVE  IS  (RD  -  R)/(V*DV> 

C 

OV  =  0.001 

C  SETTING  ALL  VD  =  V 
on  550  K  =  l.NVAR 
550  VO ( K  )  *  V ( K ) 

DO  560  J  a  1,  NVAR 
C  ADDING  DELTA  TO  VDIJ) 

VDIJ)  =  (1.  *  DV)*VIJ) 

CALL  EONSINVAR,  VD.  RD) 

DO  55B  I  ■  1,  NVAR 

POII.JJ  a  (ROII  1  -  RII  ) )/( VI J)*DV) 

558  CONTINUE 

C  RETURNING  VDIJ)  TO  VIJ)  VALUE 

VO(J)  »  VIJ) 

560  CONTINUE 
RETURN 
END 
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0001 


0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 


0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033- 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 


SUBROUTINE  GAUSSY.  (A.  ft,  x,  N) 

SOLUTION  OP  SIMULTANEOUS  EOUATI ONS  BY  GAUSS  ELIMINATION 

DIMENSION  A(81.81),x(Sl).ft(81),BC(81) 
beginning  0F  ELIMINATION  process 

00  28  K  =  l,N 

MOVING  LARGEST  COEFFICIENT  INTO  DIAGONAL  POSITION 
AnAX  s  On 
DO  4  I  a  K,N 

IFIABSIAII ,KI )  -  ABSIAMAXIJ  4.  4,  7 
2  AMAX  =  A(I,K) 

IMAX  «  I 
4  CONTINUE 

TESTING  FOR  INDEPENDENCE  OF  EOUATI ONS 
IFIABSI AMAX )  -  0.1E-15)  10.  10,  14 
10  WRITE  (6,12) 

12  FORMAT  CO  EOUATIONS  ARE  NOT  INDEPENDENT') 

Kfc TURN 

EXCHANGING  ROW  IMAX  ANO  ROW  K 
14  BTEMP  =  BIK) 

8(K)  =  B(IMAX) 

B ( IMAX  )  =  BTEMP 
DO  18  J  »  K.N 
ATEMP  a  A(  K, J ) 

A(K,J)  x  All  MAX.  J) 

18  A ( IMAX , J  )  =  ATEMP 

kSplBu?=CTJN*G  iA,,*k,/mk*m  t,mes  term  i"  fi«st  eo  from  OTHERS 

IFIK  -  N)  22,  28,  28 
22  00  24  I  a  KPLUS.N 

B(I)  =  B(I)  -  B(K)*A(I ,K)/A(K,K) 

ACON  =  A(I ,K ) 

DO  24  J  =  K.N 

24  A(I,J)  s  A ( I , J )  -  A(K,J)*ACON/A(K,K) 

28  CONTINUE 

BACK  SUBSTITUTION 
L  =  N 

32  SUM  =  0.0 

IF(L  -  N)  34,  38,  38 
34  LPLUS  =  L  +  l 

DO  36  J  =  LPLUS,  N 
36  SUM  =  SUM  +  A(L.J)*X(J) 

38  CONTINUE 

X(L)  =  CB(L)  -  SUM)/A(L.L) 

IF  (L  -  1 )  42,  42,  40 
40  L  *  L  -  1 
GO  TO  32 
42  RETURN 
END 
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0001 

0002 

0003 

0004 

000$ 

0006 

0007 

0008 

0009 

0010 

0011 

0012 

0013 

0014 

001$ 

0016 

0017 

0018 

0019 

0020 

0021 


0022 

0023 

0024 

002$ 

0026 

0027 
0028 
0029 
00  30 
0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 


SUBROUTINE  EONS (NvAR , V ,R  ) 

DIMENSION  R (8  1 ) » V ( 8 1 ) 

COMMON  A.Z.0.ALF2,ALF4.ALE6.ALF8.ALF10,BtT3.BET5.BET7.8ET9.BETll 

COMMON  H0Bl*H0B2»H063tH084,H085»Ul»U2tU3*U4*U5tAl . A3.A5.A7 »A9 

COMMON  AT.ADUCT.ALF1 

COMMON  ALTn.GRBE. ALTE 

Rl  1  )=( ALTO/I  I ALTE )* (GRBE  ) )  )-V  (  1) 

R (  2)*VI  1 )-V(30)-V(50)-V(70) 

RI  3)*l.-(.66*( ( (.028 32/. 45 36 )*(  V(  3) /Z ) )**  »28)*( 1 1 ( ( .55)*!  I  A/. 304 
18)**. 17))-. 17  )*( (68.94757/Vt 10) )  **. 22 ) )♦. 17  ) ) -V (  4) 

R<  4)*(  ( (  1020.*  ( ( 1-V  I  4)  )♦(. 075*  (VI  10/68.947 57  )) )  )*1535.  )/l  .8)-VI 
112) 

RI  5)*(<V(  2))/<<3.14159)*(iD**2.)/4.)))-VI  3) 

R (  6)*( 1 •  —  V I  4)1/12. *V  ( 4 ) )— VI  5) 

RI  7)«ll.*VI  4) )/(4.f2.*V(  4)  )-V!  6) 

RI  8)*(1.+3.*V<  4 )  )/(4.  +  2.*V(  4))-V(  7) 

RI  9  )*V(  5)*17.032+V<  6 >*28.016*V(  7)*2.016-VI  8) 

RI10)*VI 11)*1VI14)**2)/(2.*V( 15) )-V(12) 

R(ll  )=V(  9 )♦( V ( 1 3  )*V I 14)**2)/(2.E05)-V( 10) 

RI 12  )*( V( 16)*(83 14. 25/VI  8 ) )* V ( 1 1 ) ) **. 5- V( 1 4) 
Rll3)*(V(13)M83i4.25/V(  8  )  )*V(  11  )  ) / (  10.** 5) -VI  9) 

R(14)*V< 13)*V( 14)*A1-V(  2) 

RI 15 )* I (27549. 97-M  25.627418  )*( ( V( 1 1 )♦ V( 23 ) )/ 2) +  . 9900599E-02*!  I  V(  1 
11)*V<23) )/2)**2-(.668603E-05)*!(V(ll)+V(23) )/ 2) **3+ (( 55,5/1  (  V< 111* 
2V!23))/2))**2)*II IVI  9)+V(21))/2)/137.B9514)  )*VI  5)* 

31 28882. 15-( 1. 570255 )*(( VI 11)+V( 23 >)/2)+( 0.80751 2E-02)* I (V( 11)+V(23 
4 ) )/2 )**2-(  .287064E-05)*( ( V ( 1 1 ) > VC  23 ) ) /2 )** 3+ ( (  55 . 5/ 1 1 V( 11 )  +  V( 23 ) )/ 
52 ) )**2 )* I  1 1 V (  9)+V(21 ))/2)/137.89514))*Vl  6  )♦ 

6(29087.17-1 1.914598)*!  f'/(  1) >+V( 23  )  ) /2 )♦  (  0.  4001 16E-02 ) * (  (V(  11  )+V(  23 
7 ) ) / 2  )  **  2  — ( .867854E-06 )*( ( V ( 1 1 )+V(23)!/2)**3*( (55. 5/ ( ( VI 11 )  +  V(  23))/ 
82) )**2 )* ( I (V(  9)*V(21  )>/2)/46. 965050) i*V(  7))/V(  8  )-V ( 1 5  > 
RI16)=(V(15)/»V(15)-(8314.25/V(  8))))-V(16) 

R(17)=TAM( ALF2)*TAN(Vt2R ) )-TAM(  V( 27 ) )-T AN ( BET3 ) 

RI  18  )  =  BFT3-V ( 27  ) 

R(19)=2.*V(20)*U1*TAN(BET3)/10.**4-V( 19 )*V ( 1 5 )*V ( 1 2 )* ( 1 .-(  VI 22 )/V( 
110  ))**(( V(16)-l.)/V( 16)) )/10.**4 

R I  20  )*V I  2)*V ( 19 )*V( 15 )*V( 12)*( 1. -( V ( 22)/V(10) )**( (VI 16)-1. )/V( 16) 
1))/I 10.**6)-V(30) 

RI21 )=V( 12 )*( l.-VI 19)*( l-( V( 22)/V( 10) )**( ( V(  161-1.  )/V( 16)) ) )-V(24) 
RI 22  )=V( 14)*C0S(ALF2)-V(20) 

R(23)  =  .025*(l.  +  ( ( V( 27 )♦ BE T3  > /l. 57 )**2)*( l.*( 3.2/HOB1 ) )-V(17) 

R(24)  =  .025*(  l.  +  < (ALF1  +ALP2)/1.57)**2)*( 1.* (3.2 /H0B1 ) )— V( 18) 

R(25 )  =  l./( !.♦< <  V ( 18 )« I V ( 14) /U1 )**2+V( 17)*( V( 20 )/ (U1*C0S (8ET3 ) )  )**2 
l)/( I2.*V( 20) /Ul  )* I T  AN ( ALE  2 )  +  T AN( V( 28 ) ) )) ) )-V(  19  ) 

RI26)*V(21)*IV(25  )*V( 26  )*• 2 ) /( 2 . 605 ) -V( 22  ) 

RI27)=V(23)*(V(26)**2)/( 2.*V( 15) )-V(P4) 

R(28)=(V(25)*(8314,25/V(  8 ) )*V ( 23 ) ) / ( 10. ** 5) -V ( 21 ) 
R(29)=V(25)*V(26)*V(29)-V(  2) 

R(3O)=V(20)/COS(V(28 ))-V(26) 

R(31)*V(32)*(V(34 )**2 ) /( 2. * V( 35 ) ) -VI  24) 
R(32)*V(31)*(V(33)*V(34)**2)/(2.E05)-V(22) 
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0039 

0040 

0041 


0042 

0043 

0044 

0045 

0046 

0047 

0046 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 

0061 


0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 


R(  33)*( VI 33  )*(83l4'.25/VI  8)  )*VI32)  )/l 10.**5)-V(  31) 

R(34)*V(33 )*VI34)*A3-VI  2) 

RI35)*( (2754Q.97+I25.627418 )* 1 1 V  I  32 1+VI43I ) /2 )♦. 99005996-02*1 1  VI 3 
12I+VI43) )/2 1**2-1.6686036-05)* I  I  VI 32 ) +VI43 ) ) /2 )**3+l ( 55.5/11 VI 32)  + 
2VI43D/2)  )**2)*H(VI  31 1 +V 1 4 1  •)  )  /2  )  /  1 37.39514  ))  *V  I  5)  + 

31 28882. 15-1 1.570255  I* 1 1 VI 32) +VI 43) 1/21  +  1 0.  80  75  1  26-  02)* 1  1 VI 32J+VI  43 
411/21**2-1.2870646-05 1*11 VI 32 1+V 143)) /2)**3+l 1 55.5/ 1  I  V 1 32 )+VI 43 1 )  / 
52))**2)*ll  I VI 31 1+VI41 ))/2i/ 137.89514) l*VI  61  + 

6(29087.17-1 1. 914593 l*IIV( 32)  +  V(4  3) 1/21+10. 4  001166  -02)* I IVI  32I+VI43 
7  I )/2 )**?-!. 86 98546-06  )* ( I V<  32 )♦ VI 43 ) ) /? )** 3+ 1 1  55 . 5/ 1 1  VI 32 )♦ VI 43 ) ) / 
82  )  )**2  )*IHVI31)+V  141)1/2)  /45 .965050) )  *VI  7  ) ) /V  I  8I-VI35) 

R (36 >*(V(35)/(V135  1-18314. 25/VI  8)  )) )-VI36) 

RI37)*(U2/VI40)  )-<TANIAlF4)-TANIVI48))  ) 

R (38 ) *V (47 I-VI48) 

R(39)*V(40) *U2* (TANI ALF4 l+TAN ( VI 48 ) ) ) /10.**4-VI 39)*VI35)*V(  24)*(  1. 
1-(V(42)/V(22) )**I(VI36)-1.)/V(36) ))/10.**4 
R(40)*V(  2 ) *V( 39 )* VI 35 )*V(24 ) *( l.-( VI 42J/VI  22 ) ) **( I VI 361-1. )/V(  36) 

1 1 ) /( 10. **6 I-VI50) 

R(41)«VI24)*(1.-VI39 I* I  1-1 V 142) /VI 22) ) ** 1 1 V  1 36) -1 . ) /VI 36) ) ) ) — VI  44) 
R142)*V(34 )*C0SI ALF4I-VI40) 

R(43)*.025*l l.+  l IV147)  +  BET5)/1.57)**2)*I 1.+ I3.2/HCB2) I -V I 37» 

R (44)*. 025*11.  + 1 (VI 28 1  +  ALF4) /1 . 57 )** 2  I* I l. +(  3.  2/H0B2 ) )-V(33 ) 

RI45 )*!./( l.+l (V(36)*(VI34)/U2)**2+VI37)*(VI40) /(U2*COSI?ET5) ) )**2 
1 )/((2.*V(40)/U2)*(TANI6LF4)+TAN| V(4d ) ) ) ) ) J-VI39) 
R(46)*V(41)+(V(45)*V 146)** 21/12. E05I-V (42) 

R (47 ) =V(43 )  +  ( V(46  )**?)/( 2. *V 135 ) )-VI  44) 

R(48)*(V(45)*(8314.25/Vl  8  )  >* V (43 ) ) /( 10. **5  )- VI 41 ) 

R  1 49  )  =V  1 45  )  *V  1 4  6  )  *V  1 4  9  )  -V I  2) 

R(50)=V(40) /COS (V (48 ) ) -V (46 ) 

RI51 )*V(52 )♦ (V (54 )**2/( 2.*V I  55 ) ) )-V 1 44 ) 

R(S2)*V(51 )+( V(53)*V(54)**2)/<2.E05)-V(42) 

R 1 53)= I VI 53  I  * (63 14.25/V  I  8>»*V(52))/( )0.**5 )-  V(51) 

R 1 54 ) *V (53 )*V (54  ) *A5-V I  2) 

R 1 55)* 1 1 27540.97+ 1 25 .627418  )* I ( V ( 52 )  +  V 1 63 ) ) /2 )♦ .99O05O9E-02*! ! VI 5 

12) +VI63))/2)**2-( .6686036-05 )* ( I V I  52 ) +VI 46 ) ) /2  )**3+ 1  1  55 .5/1  (  VI  52)  + 
2V(63))/2))**2)*((lV(5l)+V(6l))/2)/l37.895l4))*V(  5)+ 

3 (28882.15-1 1.570255)* 1  1 V (52) +VI 63 ) )/2)  +  (  0.80 7512E-02)*! (V(  521+VI  63 
4) )/2)**2-( . 2 87 064 E- 05 )*( (V(52 )+V(63 ) )/2)**3  +  (  1 55. 5/ 1  1 V 1 52 ) +VI 63  ) ) / 
52))**2)*i; (V I  5 1 ) +V (61 ) )/2)/l 37.89514 ))* VI  6)  + 

6(29087.17-1 1 .9 14598 )*<(V(52)+V(63>>/2t+I 0.400 11 66-02)* (I VI52J+VI63 
7 ) )/2 ) **2-( .86  9  8  5  4  6-  06 )* ( I VI 52 >  +  V<  63 ) ) /2 )**3+ ( (  55 . 5/ 1  1 V<  52  >+V(  63 ) )/ 
82))**2)*(((V(5l)+V(61))/2) /45. 96 5050) )*V I  7  ) ) /V 1  ft  1— VI 55 ) 
R(56)*(V(55)/(V( 55  )-( 8 314. 25 /VI  8  ) ) ) J-VI56) 

RI57 )  =  (U3/V 160 ) >-<TAN( Al*6)-T6N< V<68  >) ) 

R  158  )=V (67  I-VI68) 

RI59)*V(60)*U3*ITAN(ALF6)  +  TAN(VI68 ) ) )/lO.**4- VI 59)*V(55)*V(  44)*(  1 . 
1-IV(62)/V(42) )**l (VI56I-1. ) / VI 56 ) ) ) / 10, **4 
R 1 60)=V I  2)*V(59)*V(55)*V144)*(l.-( VI 62) /VI  42 ))**(( VI 56) -l.) /VI  56) 
1  ))/( 10.**6)-VI70) 

R(61)*VI44)*H.-VI59)*ll-(VI62)/V(42))**UV(5fc)-l.)/V(56))))-V(  64) 
R(62)=VI54 )*COS( ALF6)-V160i 

RI63  )*.025*( 1.  + 1 IVI67  J+86T7 )/l .57 )**2)*(  1.+  1 3.2/H0B3 ) )-V(  57) 
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0070  R(64I«.025*(1.+((V(481+ALF6>/1.57>**2>*( 1.+ ( 3.2 /H063) ) -VC 58) 

0071  R(65)*l./( l.+((V(58)*(V(54)/U3)**2+V(57)*( V(60)/(U3*C0$(8ET7))>**2 

l)/( (2.*V(60)/U2 )*(TAN( ALF6)+TAN( V(68 )) ) I ) )-V< 59 > 

0072  R (66 )*V ( 61  »-*■(  V  ( 65  )«V(66)**2  )  /(  2 .605)  -V«  62  I 

0073  R(67»«V(63)+(V(66)**2)/(2.*V( 55)  I-VI64) 

0074  R(68 l*( V(65  )* (8314.25 /V (  8 ) >*V( 63 ) ) /( 10.**5)-V(6l > 

0075  R(69)«V(65!*V(66)*V(69)-V(  2) 

0076  R(70)*V(60)/C0S(V(66 ) )-V(66) 

C 

0077  R(71)XV(72)+(V(74)**2/(2.*V(79)))-V(64> 

0078  R( 72 )«V( 7 1 )  +  (VI 73 )*V (74 )**2 ) / ( 2.E05 )-V( 62  ) 

0079  R(73)*(V(73)*(8314.25/V(  8 ) >*V( 72) )/( 10.**5)-V(71 ) 

0080  R(74)«V(73)*V(74)*A0UCT-V(  2) 

0081  R(751XV(75)+(V(77 )*V( 78 )**2 > /( 2. E05)-V( 62 ) 

0082  R(76 )* ( V (  80)*(8314.25/V(  8 ) )*V( 76 )  )**. 5-V( 78 > 

0083  R(77)«(V(73)*(8314.25/V(  8 ) )*V( 76 > ) /( 10.**5> -V( 75 > 

0084  R(78)-V(77)*V(78)*V(  81)-V(  2) 

0085  R (79 )x ( ( 27549 .97* (25.627418  )* ( ( V( 7 2 )♦ V( 76 ) ) /2) ♦. 9900599E-C2*(  (  V(  7 

12)+V(76) )/2 )**2-(.668603E-05>*( ( V(72)+V(  76) >/2)**3+( ( 55 - 5/(  (  V( 72)  + 
2V( 76 ) )/2 ) )**2 )*(( (V(7i)+V(75))/2) /137. 89514)  )*V[  5  )+ 

3( 28882. 15 -( 1. 570255 )*((V(72)+V(76))/2)+( 0.807 51 2E-02)*((V( 72 )+V( 76 
4  )  )/2  )**2-( .287 064E-05  )* ( ( V( 72 )  +  V( 76) ) /2>**3+ U  55 . 5/ ( ( V( 72>+V(  76)  )/ 
52))**2)*( ((V(71)+V(75))/2)/137.89514))*v(  6 )+ 

6 (29087. 17-( 1.914598 )* ( < V ( 72 )+V ( 76 ) ) /2 )♦ ( 0. 4001 16E-02 ) * ( ( V(  72)+V(  76 
7) )/2 )**2-( .869854E-06 )* I ( V (72 ) +V (76  ) ) /2)**3+( ( 55.5/( ( V( 72  )+V( 76))/ 
82))**2)*(((V(71)+V(75))/2)/45.965050))*V(  7))/V(  8)-Vi 79) 

0086  R(80)«(V(79)/(V(79)-(8314.25/V(  8))))-V(80) 

0087  R(81)*V(76)+(V(78 )**2/( 2.*V( 79 ) ) )-V( 64) 

0088  RETURN 

0089  ENO 
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NUMBER  OF  VARIABLES  *  81 

MAXIMUM  FRACTION  CHANGE  FOR  CONVERGENCE 


O.OIOO 


VARIABLE  NUMBER  ANO  ITS  TRIAL  VALUE 


VI 

11 

9 

ROW 

9 

2.33918 

VI 

21 

m 

MFLO 

S 

2.40615 

VI 

3) 

9 

MFLX 

9 

47.48601 

VI 

4) 

9 

PHI 

9 

0.66223 

VI 

5) 

• 

XNH3 

9 

0.12687 

VI 

6) 

9 

XN2 

S 

0.31219 

VI 

71 

m 

XH2 

* 

0.56094 

VI 

81 

a 

MWT 

9 

12.03800 

VI 

91 

c 

PI 

9 

74.53757 

VI 

10) 

* 

P01 

= 

123.86876 

VI 

11) 

a 

Tl 

9 

964.45508 

VI 

12) 

s 

TOl 

9 

1120.53223 

VI 

13) 

9 

RH01 

9 

11.18985 

VI 

14) 

s 

VI 

9 

938.99512 

VI 

1$) 

9 

CPI 

9 

2824.60083 

VI 

16) 

9 

KG1 

9 

1.32366 

VI 

17) 

m 

S-Rl 

9 

0.24222 

VI 

18) 

9 

S-Nl 

9 

0.14909 

VI 

19) 

9 

EFF1 

9 

0.82226 

VI 

20) 

9 

VAX1 

S 

260.21509 

VI 

21) 

9 

P2 

9 

54.43111 

VI 

22) 

9 

P03 

9 

57.29221 

VI 

23) 

9 

T  2 

9 

950.00342 

VI 

24) 

9 

T03 

9 

962.21338 

VI 

251 

9 

RH02 

9 

8.29571 

VI 

26) 

9 

V2 

9 

262.63550 

VI 

27) 

9 

BET  2 

9 

1.03000 

VI 

28) 

= 

ALF3 

a 

-0.13588 

VI 

29) 

a 

EFA2 

a 

0.00110 

VI 

30) 

9 

P0W1 

s 

1.07600 

VI 

31) 

a 

F3 

9 

44.80257 

VI 

32) 

a 

T3 

9 

899.77344 

VI 

33) 

9 

RH03 

a 

7.20943 

VI 

34) 

9 

V3 

9 

588.62598 

VI 

35) 

9 

CP2 

9 

2774.49243 

VI 

36) 

9 

KC2 

9 

1.33144 

VI 

37) 

S 

S-RZ 

9 

0.16169 

VI 

38) 

a 

S-N2 

9 

0.13204 

VI 

39) 

9 

EFF2 

9 

0.85734 

VI 

40) 

9 

VAX2 

9 

198.01688 

VI 

41) 

a 

P4 

9 

31.11732 

VI 

42) 

= 

P05 

= 

32.20224 

VI 

43) 

= 

T4 

9 

844.66382 

VI 

44) 

= 

T05 

a 

651.99463 

VI 

45) 

s 

RH04 

9 

5.33396 

VI 

46) 

a 

V4 

9 

201.69043 

VI 

47) 

s 

RET4 

9 

0.19115 

VI 

48) 

= 

ALF5 

8 

•  0.19115 

VI 

49) 

9 

EFA4 

9 

0.00224 

VI 

50) 

9 

P0H2 

9 

0.73580 
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VI 

51) 

P5 

■ 

26.31075 

VI 

52) 

T5 

■ 

606.62500 

VI 

53) 

RH05 

m 

4.72272 

VI 

54) 

V5 

- 

499.49512 

VI 

55) 

CP3 

m 

2749.58350 

VI 

56) 

KG3 

z 

1.33545 

VI 

57) 

S-R3 

x 

0.12246 

VI 

58) 

S-N3 

z 

0.16169 

VI 

59) 

EFF3 

s 

0.86079 

VI 

60) 

X 

VAX3 

z 

166.03278 

VI 

61) 

P  6 

r 

19.80089 

VI 

62) 

P07 

■ 

20.36765 

VI 

63) 

T  6 

■ 

766.76929 

VI 

64) 

T07 

= 

772.28223 

VI 

65) 

RH06 

z 

3.73896 

VI 

66) 

V6 

• 

174.11455 

VI 

67) 

BET6 

z 

-0.26509 

VI 

68) 

ALF7 

■ 

-0.26509 

VI 

69) 

EFA6 

■ 

0.00370 

VI 

70) 

POW3 

z 

0.52738 

VI 

71) 

PEXO 

Z 

20.25191 

VI 

72) 

Tt'XO 

■ 

771.14722 

VI 

73) 

RHOX 

■ 

3.80241 

VI 

74) 

VEXD 

■ 

78.02730 

VI 

75) 

PT.iT 

z 

17.28519 

VI 

76) 

.  i 

■ 

658.18262 

VI 

77) 

P  'OT 

s 

1.00702 

VI 

78) 

VTHT 

* 

782.42798 

VI 

79) 

S 

CPEX 

m 

2682.71362 

VI 

80) 

KG6X 

u 

1.34671 

VI 

81) 

ATHT 

9 

0.00305 
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* 


variable  number  and  its  final  value 

these  values  are  for  an  alternator  output  of  2.00  MW 


VARIABLE 


VI  1)  «  P.o W  « 

VI  21  *  MFLO  * 

VI  31.  *  HrLX  = 

VI  41  «  PHI  S 

VI  51  «  X VH3  * 

VI  61  «  XN2  « 

VI  7)  «  XH2  = 

VI  8)  =  MWT  = 

VI  9)  «  Pi  = 

VI 101  «  POI  = 
VI  11)  r  T1  * 
VI 12)  «  T01  * 

VI  13)  «  RH01  = 
VI 14)  =  vi  = 
VI 15)  «  CPI  « 
VI16)  «  KG1  * 
VI 17)  *  S-Rl  * 
VI 18)  =  s-Nl  = 
VI 19)  *  Effi  , 
VI 20)  «  VAXl  * 
VI 2 1)  «  P2  * 
VI 22)  =  P03  * 

VI 23)  =  T2  * 

V 1 24 )  *  703  = 

VI 25)  *  RH02  a 
VI 26 )  »  V2  a 
VI 27)  *  8FT3 
VI 28 )  «  ALF3  = 


VI 29 )  a  EFA2  = 
VI 30)  *  P0W1  = 


VALUE  I  SI  ) 

2.33918  MW 
2.40616  KG /SEC 
47.48625  KG/SEC-SO.M 
0.66223 
0.12687 
0.31219 
0.56094 

12.03801  KG/MOLE 

74.53792  BARS 
123.86929  BARS 
964.45630  DEG  K 
1120.53320  OEG  \ 
11.1899]  KG/U.M 
938.99487  M/SEC 
2824.60352  J/KG-DEG  K 
1.32366 
0.24222 
0.14909 
0.82226 

260.21484  M/SEC 
54.43143  BARS 
57.29251  8ARS 
950.00439  OEG  K 
962.21436  OEG  K 
8.29575  KG/CU.M 
262.63550  M/SEC 
1.03000  RAOIANS 
-0.13588  RAOI ANS 
0.00110  SO.M 
1.07600  MW 


VI  31 ) 

* 

P3 

S 

VI  32) 

* 

T3 

a 

VI  33  ) 

* 

RH03 

s 

VI  34) 

* 

V3 

- 

VI35) 

9 

CP2 

= 

VI  36) 

* 

KG2 

r 

VI37) 

a 

S-R2 

- 

VI  38) 

= 

S-N2 

= 

VI 39 ) 

s 

FFF2 

r 

V 1  40) 

* 

VAX  2 

- 

V 1 4 1 ) 

* 

P4 

- 

VI  42) 

* 

P05 

- 

VI 43 ) 

= 

T  4 

- 

V 1  44 ) 

* 

T05 

r 

VI 45  ) 

* 

RHQ4 

r 

V 1  46) 

a 

V4 

_ 

VI47 ) 

9 

BET4 

VI  48) 

9 

ALF5 

a 

VI49) 

9 

EF  A4 

a 

44.80284  BARS 


899.77466 

7.20948 

588.62476 

2774.49609 

1.33144 

0.16169 

0.13204 

0.85734 

198.01649 

31.11761 

32.20251 

844.66528 

851.99609 

5.33401 

201.68987 

0.19115 

0.19115 

0.00224 


DEG  K 

KG/CU.M 

M/SEC 

J/KG-OEG  K 


M/SEC 

BARS 

BARS 

OEG  K 

OEG  K 

KG/CU.M 

M/SEC 

RAOIANS 

RAOIANS 

SO.M 


VALUE  I  BRIT ISH) 
3136.89282  HP 

5.30467  L8M/SEC 

9.72593  LBM/SEC-SO. FT 

0.66223 

0.12687 

0.31219 

0. 56094 

12.03801  LBM/MOLE 

1081.08130  PSIA 
1796.57251  PSIA 
1736.02051  OEG  R 
2016.95874  OEG  R 
0.69850  L8M/CU.FT 
3080.69189  FT/SEC 

0.67510  BTU/LBH-OEG  R 
1. 32366 
0.24222 
0. 14909 
0.82226 

853.72314  FT/SEC 
789.46118  PSIA 
830.95776  PSIA 
171D. 00708  OEG  R 
1731.98486  DEG  R 

0.51784  L8M/CU.FT 
861.66504  FT/SEC 
59.D1897  DEGREES 
-7.78581  DEGREES 
0.D1189  SO. FT 
1442.94482  HP 

649.81030  PSIA 
1619.59351  DEG  R 
0.45003  L8M/CU. FT 
1931.18359  FT/SEC 

0.66312  BTU/LBM-DEG  R 

1.33144 

0.16169 

0.13204 

0.85734 

649.66040  FT/SEC 
951. 32275  PSIA 
467.05811  PSIA 
1520.39673  DEG  R 
1533.59229  OEG  R 
0.33296  LBM/CU. FT 
661.71216  FT/SEC 
10.95272  DEGREES 
10.95272  DEGREES 
0.02407  SO. FT 
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VC  50) 

a 

POM  2 

m 

0.73580 

MM 

986.73047 

HP 

VC51) 

P5 

a 

26.31107 

BARS 

381.60962 

PSIA 

VC  52 » 

T5 

a 

806.62720  DEG  K 

1451.92822 

DEG  R 

VC  53) 

RH05 

a 

4.72277 

KG /CU.M 

0.29481 

L8M/CU.FT 

VC  54) 

V5 

m 

499.49121 

M/SEC 

1638.750^3 

FT/SEC 

VC55) 

CP3 

a 

2749.58594 

J/KG-OEG  K 

0.65717 

BTU/LBM-OEG 

R 

VC56) 

KG3 

a 

1.33545 

1.33545 

VC  57) 

S-R3 

a 

0.12246 

0.12246 

VC58) 

S-N3 

a 

0.16169 

0.16169 

VC  59) 

EFF3 

a 

0.86079 

0.86079 

VC  60 ) 

VAX3 

a 

168.03)48 

M/SEC 

551.28442 

FT/SEC 

VC  6 1 ) 

P6 

a 

19.80138 

BARS 

287.19458 

PSIA 

VC  62) 

P07 

a 

20.36813 

BARS 

295.41479 

PSIA 

VC63) 

T6 

a 

766.77124 

CEG  K 

1380.18750 

DEG  R 

VC64) 

s 

T  07 

s 

772.28394 

OEG  K 

1390.11035 

DEG  R 

VC  65) 

RH06 

a 

3.73905 

KG/CU.M 

0.23340 

L8M/CU.FT 

VC  66) 

V6 

a 

174.11438 

M/SEC 

571.24121 

FT/SEC 

VC  67) 

BET6 

a 

-0.26511 

RAOIANS 

-15.19081 

DEGREES 

VC68) 

ALF7 

s 

-0.26511 

RADIANS 

-15.19081 

DEGREES 

VC  69 ) 

EFA6 

a 

0.00370 

SO.M 

0.03978 

SO. FT 

VC  70 ) 

s 

P0W3 

a 

0.52737 

MW 

707.21655 

HP 

VC71) 

s 

PEXO 

a 

20.25240 

BARS 

293.73608 

PSIA 

VC  72) 

* 

TEXO 

a 

771.14941 

DEG  K 

1388.06812 

DEG  R 

VC73 ) 

■ 

RHOX 

a 

3.80251 

KG/CU.M 

0.23736 

L8M/CU.  FT 

VC  74) 

a 

VEXO 

a 

78.02518 

M/SEC 

255.98813 

FT/SEC 

VC75) 

8 

PTHT 

a 

17.28563 

BARS 

250.70692 

PSIA 

VC  76) 

a 

TTHT 

.a 

658.18384 

DEG  K 

1184.73022 

DEG  R 

VC  77 ) 

* 

RHOT 

a 

1.00704 

KG/CU.M 

0.06286 

IBM/CU.  FT 

VC  78) 

* 

VTHT 

a 

782.42920 

M/SEC 

2567.02515 

FT/SEC 

VC  79 ) 

s 

CPEX 

a 

2682.71533 

J /KG -DEG  K 

0.  £-4118 

8TU/L8M-DEG 

R 

VC  80 

= 

KGEX 

8 

1.34671 

1.34671 

VC  811 

* 

ATKT 

a 

0.00305 

SO.M 

0.03287 

SO. FT 

NUH06A  OF  VARIA8LES  =  81 

kaximu*  fraction  CHANGE  FOR  CONVERGENCE 


0.0100 


m 


VARIABLE  NUMBER 


VI 

1) 

* 

POU  = 

VI 

2) 

a 

MFLO  = 

VI 

3) 

s 

mflx  = 

V( 

4) 

- 

PHI  - 

VI 

5) 

s 

XNH3  « 

VI 

6) 

= 

XN2  = 

V( 

7) 

a 

XH2  = 

V( 

8) 

« 

MWT  = 

V( 

9) 

a 

PI  = 

VI 

10) 

a 

P01  = 

VI 

11) 

m 

T 1  = 

VI 

12) 

a 

T01  = 

V( 

13) 

a 

RH01  = 

VI 

14) 

3 

VI  = 

V( 

15) 

s 

CPI  = 

VI 

16) 

= 

KG1  = 

VI 

17) 

3 

S-Rl  = 

VI 

18) 

S 

S-Nl  = 

VI 

19) 

a 

EFF1  = 

V( 

20) 

a 

VAX  1  = 

VI 

21) 

a 

P2  = 

VI 

22) 

a 

P03  = 

VI 

23) 

a 

T2  = 

V(  24)  *  T 03  r 
V(  25)  a  RH02  = 


VI 

26) 

=  V2  = 

V( 

27) 

*  BET 2  = 

VI 

28) 

*  ALF3  = 

VI 

29) 

=  EFA2  = 

VI 

30) 

=  P0W1  = 

VI 

31) 

*  P3  « 

VI 

32) 

*  T3  = 

VI 

331 

*  RH03  = 

VI 

34) 

*  V3  = 

VI 

35) 

*  CP2  = 

VI 

36) 

=  KG2  = 

VI 

37) 

=  S-R2  = 

VI 

38) 

*  S-N2  = 

VI 

39) 

=  EFF2  = 

VI 

40) 

*  VAX2  = 

VI 

41) 

=  P4  = 

VI 

42) 

=  P05  = 

VI 

43) 

»  T4  = 

VI 

44) 

*  T05  = 

VI 

45) 

=  RH04  = 

VI 

46) 

-  V4  = 

VI 

47) 

“  BET4  = 

VI 

48) 

*  ALF5  = 

VI 

49) 

■  EFA4  * 

VI 

50) 

=  P0W2  a 

ANO  ITS  TRIAL  VALUE 
2.33918 
2.40616 
47.48625 
0.66223 
0.12687 
0.31219 
0.56094 
12.03801 
74.53792 
123.86929 
964.45630 
1120.53320 
11.18991 
938.99487 
2824.60352 
1.32366 
0.24222 
0.14909 
0.82226 
260.21484 
54.43143 
57.29251 
950.00439 
962.21436 
8.29575 
262.63550 
1.03000 
-0.13588 
0.00110 
1.07600 
44.80284 
899.77466 
7.20948 
588.62476 
2774.49609 
1.33144 
0.16169 
0.13204 
0.85734 
198.01649 
31.11761 
32.20251 
844.66528 
851.99609 
5.33401 
201.68987 
0.19115 
0.19115 
0.00224 
0.73580 
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VC 

511 

X 

P5  * 

26.31107 

V( 

52) 

X 

T5  = 

806.62720 

VC 

53) 

s 

RH05  * 

4.72277 

vc 

54) 

X 

V5  = 

499.49121 

VC 

55) 

s 

CP3  = 

2749.58594 

VC 

56) 

X 

KG3  = 

1.33545 

VC 

57) 

s 

S-R3  = 

0.12246 

vc 

58) 

s 

S-N3  = 

0.16169 

V( 

59) 

S 

EFF3  = 

0.86079 

V( 

60) 

s 

VAX3  = 

168.03148 

VC 

61) 

5 

P6  = 

19.80138 

VC 

62) 

r 

P  07  = 

20.36813 

VC 

63) 

X 

T  6  = 

766.77124 

VC 

64) 

X 

T07  = 

772.28394 

VC 

65) 

s 

RH06  = 

3.73905 

VC 

66) 

r 

V6  - 

174.11438 

VC 

67) 

X 

BET6  = 

-0.26511 

VC 

68) 

X 

ALF7  = 

-0.26511 

VC 

69) 

X 

EFA6  = 

0.00370 

VC 

70) 

X 

P0H3  = 

0.52737 

VC 

71) 

X. 

PEXO  = 

20.25240 

VC 

72) 

= 

TEXO  = 

771.14941 

VC 

73) 

X 

RHOX  = 

3.80251 

VC 

74) 

X 

VFXO  » 

78.02518 

VC 

75) 

X 

PTHT  s 

17.28563 

VC 

76) 

X 

TTHT  = 

658.18384 

VC 

77) 

X 

RHOT  = 

1.00704 

VC 

78) 

X 

VTHT  = 

782.42920 

VC 

79) 

x 

CPE  X  = 

2682.71533 

VC 

80) 

X 

KGEX  * 

1.34671 

VC 

81) 

.n 

ATHT  = 

0.00305 

65 


lirinriHllllMi  ■■ 


VARIABLE  NUMBER  AND  ITS  FINAL  VALUE 


THESE  VALUES  ARE  FOR  AN  ALTERNATOR  OUTPUT  OF  1.60  NW 


VARIABLE  VALUE « SI >  VALUE (BRIT ISH) 


v<  n 

Z 

PDW 

Z 

1.87134 

MW 

2509.51270 

HP 

VI  2) 

z 

MFLO 

3 

1.92066 

KG/SEC 

4.23432 

LBM/SEC 

VI  3) 

z 

MFLX 

3 

37.90465 

KG/SEC-SO.M 

7.76346 

LBM/SEC-SO.  FT 

VI  4) 

•« 

PHI 

3 

0.67791 

0.67791 

VI  51 

s 

XNH3 

= 

0.12028 

0.12028 

VI  6) 

c 

XN2 

Z 

0.31329 

0.31329 

VI  71 

z 

XH2 

3 

0.56643 

•  0.56643 

VI  8) 

— 

MWT 

S 

11.96755 

KG /MOLE 

11.96755 

LBM/MOLE 

VI  9) 

s 

PI 

S 

58.86128 

BARS 

853.71094 

PSIA 

VI10) 

= 

PD1 

= 

97.94031 

BARS 

1420.50439 

PSIA 

VI  111 

= 

T1 

= 

941.36304 

DEG  K 

1694.45264 

DEG  R 

VI 12) 

X 

T01 

3 

1095.66870 

DEG  K 

1972.20264 

OEG  R 

VI 13) 

X 

RHO 1 

= 

9.00026 

KG/CU.M 

0.56181 

LBM/CU.FT 

VI 14) 

K 

VI 

= 

931.87866 

M/SEC 

3057.34473 

FT/SEC 

VI  15) 

= 

CPI 

3 

2813.89038 

J/KG-DEG  K 

0.67254 

BTU/LBM-OEG 

R 

VI 16) 

z 

KG1 

= 

1.32783 

1.32783 

VI  17) 

= 

S-Rl 

= 

0.24222 

0.24222 

VI 18) 

z 

S-Nl 

3 

0.14909 

0.14909 

VI  19) 

z 

EFF1 

3 

0.82337 

0.82337 

VI 20) 

s 

VAX1 

3 

258.24268 

M/SEC 

B47. 25293 

FT/SEC 

VI  21) 

= 

P2 

3 

42.71643 

BARS 

619.54956 

PSIA 

V 1  22 ) 

c 

P03 

= 

44.97217 

BARS 

652.26611 

PSIA 

VI  23) 

= 

T2 

3 

925.88037 

OEG  K. 

1666.58374 

DEG  R 

VI  24) 

s 

T  03 

3 

937.95166 

06G  K 

1688.31226 

OEG  R 

VI  25) 

= 

RHO  2 

3 

6.64084 

KG/CU.M 

0.41453 

LBM/CU.FT 

V 1  26 ) 

3 

V  2 

= 

260.64502 

M/SEC 

855.13452 

FT/SEC 

VI  27) 

= 

BET2 

S 

1.03000 

RADIANS 

59.01897 

DEGREES 

V 1 28 ) 

s 

ALF3 

3 

-0.13588 

RADIANS 

-7.78581 

DEGREES 

VI  29) 

= 

E  FA  2 

3 

0.00111 

SO.M 

0.01194 

SO. FT 

VI 30 ) 

z 

P0W1 

3 

0.85238 

MW 

1143.06470 

HP 

VI 31 ) 

z 

P3 

= 

35.00558 

BARS 

507.71313 

PSIA 

VI  32) 

z 

T3 

3 

875.31250 

DEG  K 

1575.56177 

DEG  R 

VI  33) 

z 

RH03 

= 

5.75647 

KG/CU.M 

0.35933 

LBM/CU.FT 

VI  34) 

= 

V  3 

3 

588.45093 

M/SEC 

1930.61328 

FT/SEC 

VI 35 ) 

X 

CP2 

3 

2764.04565 

J/KG-DEG  K 

0.66062 

BTU/LBM-DEG 

R 

VI 36) 

= 

KG  2 

3 

1.33573 

1.33573 

VI37) 

3 

S-R2 

3 

0.16161 

0.16161 

VI  38) 

- 

S-N2 

3 

0.13204 

0.13204 

V  1 39 ) 

= 

EFF2 

3 

0.85737 

0.85737 

VI  40) 

= 

VAX2 

3 

197.95799 

M/SEC 

649.46851 

FT/SEC 

V 1 4 1 ) 

X 

P4 

= 

24.10658 

BARS 

349.63647 

PSIA 

VI  42) 

- 

P05 

X 

24.96648 

BARS 

362.10815 

PSIA 

V  1 43 ) 

X 

T  4 

3 

820.02588 

DEG  K 

1476.04590 

DEG  R 

V I  44 ) 

3 

T05 

3 

827.37769 

DEG  K 

1489.27905 

DEG  R 

V ( 45 ) 

3 

RH04 

= 

4.23147 

KG/CU.M 

0.26414 

LBM/CU. FT 

VI  46) 

S 

V4 

3 

201.60138 

M/SEC 

661.42188 

FT/SEC 

VI 47) 

3 

BET4 

= 

0.19040 

RADI  ANS 

10.91018 

DEGREES 

VI  48) 

3 

ALF5 

3 

0.19040 

RADIANS 

10.91018 

DEGREES 

VI49) 

8 

EFA4 

3 

0.00225 

SO.M 

0.02423 

SO.  FT 
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V(  50) 

x 

POM  2 

X 

0.58701 

MW 

787.19580 

HP 

V(  5 11 

X 

P5 

- 

20.20836 

EARS 

293.09761 

PSIA 

VC52) 

X 

T5 

X 

780.69165 

OEC-  K 

1605.24614 

DEG  R 

V(53) 

X 

R«05 

= 

3.72592 

KG /CU.M 

0.23258 

LBM/CU.FT 

VC54) 

C 

V5 

X 

505.37764 

M/SEC 

1658.06250 

FT/SEC 

VC  55) 

* 

CP3 

r 

2735.35205 

J/KG-DEG  K 

0.65376 

BTU/L8M-DEG 

R 

v(56) 

s 

KG3 

= 

1.36065 

1.36065 

VC  57) 

s 

S-R3 

= 

0.12483 

0.12483 

VC58) 

r 

S-N3 

= 

0.16161 

0. 16161 

VC  59) 

X 

EFF3 

= 

0.86010 

0.86010 

VC  60 ) 

X 

VAX3 

170.01160 

M/SEC 

557.78076 

FT/SEC 

VC  6 1) 

X 

P6 

= 

16.95262 

BARS 

216.86946 

PSIA 

VI62) 

r 

P07 

= 

15.39656 

BARS 

223.30832 

PSIA 

VC63) 

X 

T6 

= 

739.58252 

DEG  K 

1331.24780 

DEG  R 

VC  66) 

= 

T  07 

= 

745.15942 

OEG  K 

1341.28638 

DEG  R 

VC65) 

s 

RH06 

= 

2.91013 

KG /CU.M 

0.18166 

U8M/CU.FT 

VC  66) 

s 

V6 

s 

174.67152 

M/SEC 

573.06934 

FT/SEC 

VC  67) 

= 

BET  6 

= 

-0.23151 

AO  IANS 

-13.26534 

OEGREES 

VC  68 ) 

X 

AIF7 

= 

-0.23151 

RADI  ANS 

-13.26534 

OEGREES 

VI  69) 

s 

EFA6 

X 

0.00378 

SO  .M 

0.04067 

SO. FT 

VC  70) 

= 

P0W3 

s 

0.43195 

MW 

579.25098 

HP 

VC  7 1 ) 

x 

PEXO 

S 

15.30184 

BARS 

221.93449 

PSIA 

VC  72) 

X 

TEXD 

= 

743.96313 

OEG  K 

1339.13306 

OEG  R 

VC73) 

M 

RHOX 

X 

2.96056 

KG/CU.M 

0.18480 

L8M/CU. FT 

VC74) 

X 

VEXO 

X 

79.99350 

M/SEC 

262.44580 

FT/SEC 

VC  75 ) 

= 

PTHT 

= 

13.03802 

8  ARS 

189.10060 

PSIA 

VC  76) 

M 

TTHT 

s 

633.89319 

DEG  K 

1141.01611 

DEG  R 

VC  77 ) 

= 

RHOT 

= 

0.79281 

KG/CU.M 

0.04949 

L8M/CU. FT 

VC  78) 

= 

VTHT 

= 

771.35156 

M/SEC 

2530.68115 

FT/SEC 

VC  79 ) 

= 

CPEX 

X 

2673.81494 

J/KG-DEG  K 

0.63906 

8TU/LBM-DEG 

R 

VI  80) 

1  X 

KGEX 

r 

1.35104 

1.35104 

VI  81) 

1  x 

ATHT 

= 

0.00314 

SO.M 

0.03381 

SO.  FT 
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NUMBER  Of  VARIABLES  >  81 

MAXIMUM  FRACTION  CHANGE  FOR  CONVERGENCE  - 


VARIABLE 

number  and 

ITS  TRIAL  VALUE 

VI 

1* 

* 

ROW 

9 

1.87134 

VI 

2) 

B 

MFLO 

9 

1.92066 

V( 

31 

m 

MFLX 

9 

37.90465 

VI 

41 

9 

PHI 

9 

0.67791 

VI 

Si 

a 

XNH3 

9 

0.12028 

V! 

6) 

9 

XN2 

9 

0.31329 

VI 

71 

a 

XH2 

9 

0.56643 

VI 

8) 

9 

MWT 

9 

11.96755 

VI 

9) 

9 

PI 

9 

58.86128 

VI 

10) 

9 

P01 

9 

97.94031 

VI 

11) 

9 

T1 

9 

941.36304 

VI 

12) 

9 

T01 

X 

1095.66870 

VI 

13) 

9 

RH01 

9 

9.00026 

VI 

14) 

9 

VI 

9 

931.87866 

VI 

15) 

9 

CPI 

9 

2813.89038 

VI 

16) 

9 

KG1 

9 

1.32783 

VI 

17) 

9 

S-Rl 

9 

0.24222 

VI 

18) 

9 

S-Nl 

9 

0.14909 

VI 

19) 

m 

EFF1 

9 

0.82337 

VI 

20) 

9 

VAX  1 

= 

258.24268 

VI 

21) 

9 

P2 

9 

42.71643 

VI 

22) 

9 

P03 

9 

44.97217 

VI 

23) 

9 

T2 

9 

925.88037 

VI 

24) 

9 

T03 

9 

937.95166 

VI 

25) 

9 

RH02 

9 

6.64084 

VI 

26) 

9 

V2 

9 

260.64502 

VI 

27) 

9 

8ET2 

9 

1.03000 

VI 

28) 

9 

ALF3 

9 

-0.13588 

VI 

29) 

9 

EFA2 

X 

0.00111 

VI 

30) 

9 

P0W1 

S 

0.85238 

VI 

31) 

9 

P3 

9 

35.00558 

VI 

32) 

9 

T3 

9 

875.31250 

VI 

33) 

9 

RH03 

S 

5.75647 

VI 

34) 

9 

V3 

B 

588.45093 

VI 

35) 

9 

CP2 

r 

2764.04565 

VI 

36) 

9 

KG2 

8 

1.33573 

VI 

37) 

9 

S-R2 

s 

0.16161 

VI 

38) 

9 

S-N2 

s 

0.13204 

VI 

39) 

9 

EFF2 

s 

0.85737 

VI 

40) 

9 

VAX2 

3 

197.95799 

VI 

41) 

9 

P4 

S 

24.10t>58 

VI 

42) 

9 

P05 

X 

24.96648 

VI 

43) 

9 

T4 

9 

820.02588 

VI 

44) 

9 

T05 

r 

827.37769 

VI 

45) 

9 

RH04 

9 

'  4.23147 

VI 

46) 

S 

V4 

= 

201.60138 

VI 

47) 

s 

8ET4 

3 

0.19040 

VI 

48) 

9 

ALF5 

X 

0.19040 

VI 

49) 

9 

EFA4 

S 

0.00225 

VI 

50) 

9 

P0H2 

X 

0.58701 

O.OIOO 
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VI 

51) 

x 

P5 

X 

VI 

52) 

x 

T5 

* 

VI 

53) 

x 

RH05 

X 

VI 

54) 

X 

V5 

X 

VI 

55) 

x 

CP3 

X 

VI 

56) 

a 

KG3 

X 

VI 

57) 

s 

S-R3 

X 

VI 

58) 

* 

S-N3 

X 

VI 

59) 

s 

EFF3 

S 

VI 

6  05 

m 

VAX3 

X 

VI 

61) 

m 

P6 

X 

VI 

62) 

x 

P07 

X 

VI 

63) 

X 

T6 

X 

VI 

64) 

m 

T07 

X 

VI 

65) 

X 

RH06 

s 

VI 

66) 

a 

V6 

r 

VI 

67) 

* 

BET6 

a 

VI 

68) 

X 

ALF7 

X 

VI 

69) 

r 

EFA6 

X 

VI 

70) 

X 

P0W3 

s 

VI 

71) 

= 

PEXD 

X 

VI 

72) 

* 

TEXO 

r 

VI 

73) 

m 

RHOX 

s 

VI 

74) 

s 

VEXO 

a 

VI 

75) 

m 

PTHT 

X 

VI 

76) 

X 

TTHT 

X 

VI 

77) 

= 

RHOT 

S 

VI 

78) 

X 

VTHT 

X 

VI 

79) 

s 

CPEX 

3 

VI 

80) 

X 

KGEX 

X 

VI 

81) 

a 

ATHT 

X 

20.20836 

780.69165 

3.72592 

505.37744 

2735.35205 

1.34045 

0.12483 

0.16161 

0.86010 

170.01160 

14.95262 

15.39656 

739.58252 

745.15942 

2.91013 

174.67152 

-0.23151 

-0.23151 

0.00378 

0.43195 

15.30184 

743.96313 

2.96056 

79.99350 

13.03802 

633.89819 

0.79281 

771.35156 

2673.81494 

1.35104 

0.00314 


variable  number  ano  its  final  value 


THESE  VALUES  ARE  FOR  AN  ALTERNATOR  OUTPUT  OF  1.20  MW 


VARIABLE 
V (  1 »  *  POW  * 

VI  2)  *  melo  * 

VI  31  =  HfLX  = 

VI  4)  *  PHI  * 

VI  51  *  XNH3  * 

VI  6)  *  XN2  * 

VI  7)  *  XH2  = 

VI  8)  «  MWT  = 


VALUE  IS!  ) 

1.40351  MW 
1.43582  KG/SEC 
28.33632  KG/SEt-SO.M 
0.69686 
0.11240 
0.31460 
0.57300 

11.88343  KG/MOLE 


VALUE (BRITISH) 
1882.13501  HP 

3.16544  LBM/SEC 

5.80372  LBM/SEC -SO. FT 

0.69686 

0.11240 

0.31460 

0.57300 

11.88343  LBM/MOLE 


VI  9) 

X 

PI 

VI 10) 

* 

P01 

Vlll) 

St 

T1 

VI 12 ) 

x 

T01 

VI 13) 

St 

RH01 

VI 14) 

= 

VI 

VI 15) 

St 

CPI 

VI 16) 

= 

KG1 

VI  17) 

St 

S-Rl 

VI 18) 

St 

S-Ml 

VI  19) 

x 

PFF1 

VI  20) 

x 

VAXl 

VI  21) 

X 

P2 

VI22) 

X 

P03 

VI  23) 

s 

T2 

V 1 24 ) 

r 

T03 

VI  25) 

X 

RH02 

VI  26) 

X 

V2 

VI  27) 

x 

BFT2 

VI 28 ) 

St 

ALF3 

VI  29) 

X 

EFA2 

VI  30 ) 

c 

P0W1 

43.50079  BARS 
72.48459  BARS 
916.79028  OEG  K 
1069.23584  OEG  K 
6.78182  KG/CU.M 
924.52661  M/SEC 
2803.45947  J/KG-DEG  K 
1.33256 
0.24222 
0.14909 
0.82451 

256.20557  M/SEC 
31.35936  BARS 
33.02396  BARS 
900.25488  OEG  K 
912.18066  OEG  ** 
4.97875  KG/CU.M 
258.58887  M/SEC 
1.03000  RADIANS 
-0.13588  R  AO  I  AN S 
0.00112  SO.M 
0.63219  MW 


630.92578 
1051.30029 
1650.22168 
1924.62354 
0.42333 
3033.22388 
0.67004 
1.33256 
0.24222 
0. 14909 
0.82451 
840.56934 
454.82910 
478.97192 
1620.45001 
1641.92432 
0.31078 
848.38867 
59.01897 
-7.78581 
0.01200 
847.77856 


PSIA 
PSIA 
OEG  R 
OEG  R 
LBM/CU.FT 
FT/SEC 

BTU/LBM-DEG  R 


FT/SEC 

PSIA 

PSIA 

OEG  R 

OEG  R 

LBM/CU.FT 

FT/SEC 

DEGREES 

DEGREES 

SO. FT 

HP 


VI 31 ) 

=• 

P3 

= 

25.57379 

BARS 

VI  32) 

C 

T3 

= 

849.32617 

OEG  K 

VI  33 ) 

X 

RH03 

= 

4.30367 

KG/CU.M 

VI  34) 

•= 

V3 

= 

588.40820 

M/SEC 

VI35) 

X 

CP2 

= 

2754.15894 

J/KG -OEG 

VI  36) 

X 

KG2 

= 

1.34054 

VI 37 ) 

X 

S-R2 

= 

0.16159 

VI  38) 

X 

S-N2 

= 

0.13204 

VI 39 ) 

= 

EFF2 

= 

0.85738 

VI  40) 

= 

VAX  2 

= 

197.94360 

M/SFC 

V 1 4 1 1 

s 

P4 

17.44447 

BARS 

VI  42) 

= 

P05 

= 

18.08258 

BAP.S 

VI43) 

X 

T  4 

= 

793.84839 

OEG  K 

VI  44) 

X 

T05 

= 

801.22510 

OEG  K 

V|«5) 

X 

RH04 

= 

3.14079 

KG/CU.M 

VI  4c  ) 

X 

V4 

X 

201.57959 

M/SEC 

VI47 1 

X 

RFT4 

X 

0.19022 

RAO I  ANS 

VI  48) 

X 

ALF5 

St 

0.19022 

RAOIANS 

V(49) 

X 

EFA4 

=• 

0.00227 

SO.M 

370.91650  PSIA 
1528.78638  OEG  R 

0.26864  LBM/CU.FT 
1930.47314  FT/SEC 

0.65826  BTU/LBM  -DEG  R 

1.34054 

0.16159 

0.13204 

0.85738 

649.42114  FT/SEC 
253.01076  PSIA 
262.26563  PSIA 
1428.92627  OEG  R 
1442.20435  OEG  R 
0.19605  L8M/CU. FT 
661.35034  FT/SEC 
10.89970  OEGREES 
10.89970  DEGREES 
0.02441  SO.FT 
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VC  5  01 

X 

POM  2 

X 

0.43877 

MW 

588.40308 

HP 

V(51! 

* 

P5 

X 

14.47696 

BARS 

209.97063 

PSU 

V(52) 

r 

T  5 

X 

753.01855 

OEG  K 

1355.43262 

OEG  R 

V(  53) 

X 

P.H05 

X 

2.74783 

KG/CU.M 

0.17153 

LBM/CU.FT 

VC54I 

r 

V5 

X 

512.28345 

M/SEC 

1680.71997 

FT/SEC 

VC  55) 

s 

CP3 

3 

2721.97559 

J/KG-OEG  r 

0.6505? 

BTU/LBM-OEG  R 

VC56) 

* 

KG3 

X 

1.34596 

1.34596 

VC  57) 

S-R3 

s 

0.12770 

0.12770 

VC  58) 

S-N3 

X 

0.16159 

0.16159 

VC  59) 

= 

EFF3 

r 

0.85918 

0.85918 

VC  60) 

VAX3 

X 

172.33487 

M/SEC 

565.40308 

FT/SEC 

VC  61) 

P6 

a 

10.49730 

BARS 

152.25055 

PSIA 

VC  62) 

P07 

s 

10.82280 

BARS 

156.97144 

PSIA 

VC  63) 

T6 

3 

710.47461 

OEG  K 

1278.85352 

DEG  R 

VC  64 ) 

3 

T07 

s 

716.13696 

OEG  K 

1289.04590 

OEG  R 

VC  65) 

RH06 

a 

2.11178 

KG /CU.M 

0.13182 

LBM/CU.FT 

VC  66) 

V6 

X 

175.57477 

M/SEC 

576.03271 

FT/SEC 

VC  67) 

BFT6 

St 

-0.19241 

RADIANS 

-11.02494 

OEGREES 

VC68I 

X 

Alr7 

X 

-0.19241 

RADIANS 

-11.02494 

OEGREES 

VC  69) 

3 

FFA6 

X 

0.00387 

SO.M 

0.04168 

SO. FT 

VC  70) 

X 

P0H3 

X 

0.33255 

MW 

445.95313 

HP 

VC  7 1 ) 

x 

PFXO 

X 

10.74988 

8  ARS 

155.91385 

PSIA 

VC  721 

X 

TEXO 

X 

714.86450 

OEG  K 

1286.75537 

OEG  R 

VC  73  1 

= 

RHOX 

X 

2.14931 

KG/CU.M 

0.13416 

LBM/CU.  FT 

VC  74) 

X 

VEXD 

X 

82.37239 

M/SEC 

270.25049 

FT/SEC 

VC75 ) 

3 

PTHT 

X 

9.14250 

BARS 

132.60077 

PSIA 

VC  76) 

3 

TTHT 

X 

607.97388 

OEG  K 

1094.35229 

DEG  R 

VC  77 ) 

3 

RHOT 

X 

0.58271 

KG/CU.M 

0.03637 

LBM/CU. FT 

VC  78) 

X 

VTHT 

X 

759.42236 

M/SEC 

2491.54321 

FT/SEC 

VC  79 ) 

X 

CPEX 

X 

2665.98145 

J/KG-DEG  K 

0.63718 

BTU/LBM-OEG  R 

VC  80) 

3 

KGEX 

X 

1.35581 

1.35581 

VC  811 

X 

AThT 

X 

0.00324 

SO.M 

0.03493 

•SO.  FT 

NUMBER  of  VARIABLES  = 
maximum  fraction  change 


81 

for  convergence  * 


VARIABLE  NUMBER  ANO 
VI  1)  a  POW  * 

VI  .  2)  *  MFLO  = 

VI  3)  ■  MFLX  a 

VI  4)  -  phi  a 

VI  51  *  XNH3  a 

VI  6)  *  XN2  * 

VI  7)  «  XH2  - 
VI  8)  a  MWT  a 
VI  91  «  Pi  « 

VI  10)  »  poi  * 

VI  11)  =  Tl  . 

VI  12)  ■  T01  a 

VI  13)  «  RH01  « 

VI  14)  a  VI  = 

VI  15)  •  CPI  = 

VI  16)  *  KGl  a 
VI  17)  .  S-R 1  = 

VI  18)  -  S-Nl  = 

VI  19)  •  EFFl  a 
VI  20)  a  VAXl  a 
VI  21)  *  P2  a 
VI  22)  a  PQ3  > 

VI  23)  a  T2  a 
VI  24)  a  T03  a 
VI  25)  a  RH02  a 
VI  26)  a  V2  a 
VI  27)  a  BET2  a 
VI  28)  *  ALF3  = 

VI  29)  a  EFA2  a 
VI  30)  a  P0W1  a 
VI  31)  «  P3  = 


VI  32)  a  T3 
VI  33)  a  RH03 
VI  34)  a  V3 
VI  35)  a  CP2 
VI  36)  a  KG2 
VI  37)  =  S-R 2 
VI  38)  a  S-N2 
VI  39)  =  EFF2 
VI  40)  a  VAX2 


VI  41) 
VI  42) 
VI  43) 


P4 

P05 

T4 


VI  44)  =  T05  a 

VI  45)  a  RH04  a 
V  I  4  6 )  a  V4  a 


VI  47)  a  8ET4  a 
VI  48)  a  ALF5  a 
VI  49)  a  EFA4  a 
VI  50)  a  P0W2  a 


ITS  trial  VALUE 
1.40351 
1.43582 
28.33632 
0.69686 
0.11240 
0.31460 
0.57300 
11.88343 
43.50079 
72.48459 
916.79028 
1069.23584 
6.78182 
924.52661 
2803.45947 
1.33256 
0.24222 
0.14909 
0.82451 
256.20557 
31.35936 
33.02396 
900.25488 
912.18066 
4.97875 
25o. 58887 
1.03000 
-0.13588 
0.00112 
0.63219 
25.57379 
849.32617 
4.30367 
588.40820 
2754.15894 
1.34054 
0.16159 
0.13204 
0.85738 
197.94360 
1 7 .4444  7 
18.08258 
793.84839 
801.22510 
3.14079 
201.57959 
0.19022 
0. 19022 
0.00227 
0.43877 


0.0100 


VC 

51) 

9 

P5 

8 

14.47696 

VC 

52) 

9 

T5 

9 

753.01855 

VC 

53) 

9 

RH  5 

9 

2.74783 

VC 

54) 

9 

V5 

9 

512.28345 

VC 

55) 

r. 

CP3 

3 

2721.97559 

VC 

56) 

9 

KG3 

* 

1.34596 

VC 

57) 

9 

S-R3 

3 

0.12770 

VC 

56) 

8 

S-N3 

X 

0.16159 

VC 

59) 

9 

EFF3 

9 

0.85918 

VC 

60) 

S 

VAX3 

8 

172.33487 

VC 

61) 

8 

P6 

= 

10.49730 

VC 

62) 

8 

P07 

= 

10.82280 

VC 

63) 

* 

T6 

8 

710.47461 

VC 

64) 

9 

T07 

8 

716.13696 

VC 

65) 

9 

RH06 

9 

2.11178 

VC 

66) 

9 

V6 

3 

175.57477 

VC 

67) 

9 

BET6 

X 

-0.19241 

VC 

68) 

9 

ALF7 

8 

-0.19241 

VC 

69) 

3 

EFA6 

= 

0.00387 

VC 

70) 

S 

P0W3 

S 

0.33255 

VC 

71) 

r 

PEXO 

8 

10.74988 

VC 

72) 

8 

TEXO 

8 

714.86450 

VC 

73) 

9 

RHOX 

8 

2.14931 

VC 

74) 

9 

VEXO 

8 

82.37239 

VC 

75) 

9 

PTHT 

8 

9.14250 

V( 

76) 

S 

TTHT 

S 

607.97388 

VC 

77) 

= 

RHOT 

8 

0.58271 

VC 

78) 

8 

VTHT 

= 

759.42236 

VC 

79) 

r 

CPEX 

X 

2665.98145 

VC 

80) 

8 

KGEX 

9 

1.35581 

VC 

81) 

8 

ATHT 

9 

0.00324 
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VARIABLE  NUMBER  AMO  ITS  FINAL  VALUE 


fc 

l 

* 

i 

t 


THESE 

VALUES 

ARF 

FOR  AN  ALTERNATOR  OUTPUT  OF 

0.80  MM 

VARIABLE 

VALUE  I  SI  ) 

VALUE  1  BRITISH) 

V(  1) 

C 

POW 

a 

0.93567 

MM 

1254.75757 

HP 

V(  2) 

* 

MFLO 

= 

0.95267 

KG/SEC 

2.10028 

LBM/SEC 

V(  3)' 

= 

MFLX 

a 

18.80127 

KG/SEC-SO.M 

3.85079 

L8M/SEC-S0. 1 

VI  41 

* 

PHI 

a 

0.72132 

0.72132 

VI  51 

a 

XNH3 

= 

0.10240 

0.10240 

V(  61 

X 

XN2 

= 

0.31627 

0.31627 

VI  71 

s 

XH2 

= 

0.58133 

0.58133 

VI  8) 

* 

MKT 

X 

11.77661 

KG/MOLE 

11.77661 

LBM/MOLE 

VI  91 

a 

PI 

X 

28.49974 

8ARS 

413.35376 

PSIA 

VI  101 

X 

P01 

= 

47.56911 

BARS 

689.93164 

PS1A 

VI  11) 

X 

T1 

= 

889.58203 

OEG  K 

1601.24683 

OEG  R 

VI 12) 

a 

T01 

a 

1040.01636 

OEG  K 

1872.02856 

DEG  R 

VI 13) 

* 

RH01 

a 

4.53787 

KG /CU.M 

0.28326 

LBM/CU.FT 

VI 14) 

a 

.  VI 

a 

916.76270 

M/SEC 

3007.75171 

FT/SEC 

VI  15) 

s 

CPI 

= 

2793.42432 

J/KG-OEG  K 

0.66764 

BTU/LBM-DEG 

VI 16) 

X 

KG1 

= 

1.33821 

1.33821 

VI 17) 

X 

S-Rl 

a 

0.24222 

0.24222 

VI 18) 

X 

S-Nl 

s 

0.14909 

0. 14909 

VI 19) 

X 

EFF1 

s 

0.82573 

0.82573 

VI  20) 

X 

VA  XI 

= 

254.05411 

M/SEC 

833.51074 

FT/SEC 

VI  21) 

X 

P2 

= 

20.39499 

8ARS 

295.80420 

PSIA 

VI  22 ) 

X 

P03 

a 

21.48415 

BARS 

311.60107 

PSIA 

VI  23) 

X 

T2 

= 

871.95215 

OEG  K 

1569.51318 

OEG  R 

VI  24) 

X 

T03 

a 

883.72070  OEG  K 

1590.69653 

OEG  ft 

VI  25) 

s 

RH02 

= 

3.31306 

KG /CU.M 

0.20681 

LBM/CU.FT 

VI  26) 

X 

V2 

a 

256.41748 

M/SEC 

841.26465 

FT/SEC 

VI  27) 

X 

BET2 

r 

1.03000 

RADIANS 

59.01897 

OEGREES 

VI 28 ) 

X 

ALF3 

a 

-0.13588 

RADI  ANS 

-7.78581 

DEGREES 

VI  29) 

a 

EFA2 

a 

0.00112 

SO  .M 

0.01207 

SO. FT 

VI  30) 

X 

P0W1 

a 

0.41594 

MM 

557.78052 

HP 

VI 3 1 ) 

= 

P3 

a 

16.53975 

8ARS 

239.88884 

PSIA 

VI 32) 

a 

T3 

a 

820.62891 

OEG  K 

1477.13135 

DEG  R 

VI  33 ) 

X 

RH03 

a 

2.85483 

KG/CU.M 

0.17820 

LBM/CU.  FT 

VI  34) 

X 

V3 

3 

588.54761 

M/SEC 

1930.93042 

FT/SEC 

VI 35 ) 

a 

CP2 

a 

2745.11304 

J/KG-DEG  K 

0.65610 

BTU/LBM-DEG 

VI  36) 

X 

KG2 

a 

1.34623 

1.34623 

VI  37 ) 

X 

S-R2 

a 

0.16165 

0.16165 

VI  38) 

X 

S-N2 

= 

0.13204 

0.13204 

VI 39) 

X 

FFF2 

a 

0.85735 

0.85735 

VI  40) 

= 

VAX2 

a 

197.99052 

M/SEC 

649.57520 

FT/SEC 

V 1 4  1 ) 

X 

P  4 

a 

11.15817 

BARS 

161.83556 

PSIA 

VI  42) 

X 

P05 

a 

11.57824 

BARS 

167.92827 

PSIA 

VI 43 ) 

X 

T  4 

a 

764.94360 

OEG  K 

1376.89771 

DEG  R 

VI  44) 

X 

T05 

a 

772.34985 

DEG  K 

1390.22900 

DEG  R 

VI45) 

X 

RH04 

a 

2.06615 

KG/CU.M 

0.12897 

LBM/CU.FT 

VI  46) 

s 

V4 

a 

201.65057 

M/SEC 

661.58325 

FT/SEC 

VI47 ) 

X 

BET4 

X 

0.19082 

RADIANS 

10.93385 

OEGREES 

VI  48 ) 

X 

ALF5 

a 

0.19082 

RADIANS 

10.93385 

OEGREES 

VI49) 

X 

EFA4 

a 

0.00229 

SO.M 

0.02461 

SO. FT 

74 


t 

VI 501 

* 

POM  2 

x 

0.29126 

MM 

390.58130 

HP 

i 

VI51) 

P5 

x 

9.14639 

BAR  S 

132.65717 

PSIA 

| 

VI52» 

T5 

X 

722.31396 

OEG  K 

1300.16455 

DEG  R 

VI 53) 

RH05 

X 

1.79358 

KG/CU.M 

0.11196 

18M/CU.FT 

i 

VI  54) 

V5 

X 

520.74292 

M/SEC 

1708.47412 

FT/SEC 

f 

£ 

VI 55) 

CP3 

X 

2709.78857 

J/KG-DEG  K 

0.64766 

BTU/LBM-DEG 

R 

?■ 

VI  56) 

KG3 

s 

1.35233 

1.35233 

VI  57) 

S-R3 

X 

0.13131 

0. ' 3131 

i 

VI 58 ) 

S-N3 

X 

0.16165 

0.16165 

! 

VI  59) 

EFF3 

s 

0.85791 

0.85791 

Jp 

VI60) 

VAX3 

X 

175.18066 

M/SEC 

574.73975 

FT/SEC 

| 

VI  61) 

P6 

X 

6.46630 

8ARS 

93.78574 

PSIA 

jj 

VI 62 ) 

P07 

X 

6.67799 

BARS 

96.85611 

PSIA 

y- 

VI63) 

16 

X 

678.06128 

OEG  K 

1220.50977 

OEG  R 

VI  64) 

T07 

X 

683.84497 

OEG  K 

1230.92041 

DEG  R 

f 

VI  65) 

RH06 

X 

1.35078 

KG/CU.M 

0.08432 

IBH/CU.FT 

i 

VI 66 ) 

V6 

X 

177.04234 

M/SEC 

580.84741 

FT/SEC 

f 

VI  67) 

8ET6 

X 

-0.14515 

RAO TANS 

-e.M713 

DEGREES 

t 

VI 68 ) 

ALF7 

X 

-0.14515 

RADI ANS 

-8.31713 

DEGREES 

< 

VI  69) 

FFA6 

X 

0.00398 

SO.M 

0.04288 

SO. FT 

VI 70) 

P0W3 

X 

0.22848 

MW 

306.39526 

HP 

VI  71 ) 

PE  XD 

s 

6.62783 

8ARS 

96.12863 

PSIA 

VI72) 

TEXD 

X 

682.47412 

DEG  K 

1228.45288 

OEG  R 

VI 73 ) 

RHOX 

X 

1.37557 

KG/CU.M 

0.08587 

L8M/CU. FT 

VI 74) 

= 

VEXD 

= 

85.39653 

M/SEC 

280.17236 

FT/SEC 

VI  75) 

PTHT 

X 

5.62490 

8ARS 

81.58223 

PSIA 

VI76) 

TTHT 

X 

579.20068 

DEG  K 

1042.56055 

DEG  R 

VI 77 ) 

RHOT 

X 

0.37836 

KG/CU.M 

0.02362 

LBM/CU. FT 

VI 78) 

= 

VTHT 

s 

746.10352 

M/SEC 

2447.84619 

FT/SEC 

VI79) 

X 

CPEX 

X 

2659.82422 

J /KG -OEG  K 

0.63571 

BTU/LBM-DEG 

R 

VI  80) » 

KGEX 

X 

1.36134 

1.36134 

VI  81)« 

ATHT 

X 

0.00337 

SO.M 

0.03633 

SO.  FT 

NUMBER  Of  VARIABLES  •  81 

MAXIMUM  FRACTION  CHANCE  FOR  CONVERGENCE  ■ 


VARIABLE 

NUMB 

VI 

11 

9 

ROM 

VI 

2) 

■ 

MFLO 

VI 

3) 

u 

MFLX 

VI 

41 

u 

PHI 

VI 

51 

m 

XNH3 

VI 

6) 

u 

XN2 

VI 

7) 

m 

XH2 

VI 

81 

9 

MWT 

VI 

91 

u 

PI 

VI 

101 

w 

P01 

VI 

111 

u 

T1 

VI 

12) 

9 

T01 

VI 

13) 

u 

RH01 

VI 

14) 

u 

VI 

VI 

15) 

9 

CPI 

VI 

16) 

9 

KG1 

VI 

17) 

9 

S-Rl 

VI 

18) 

9 

S-Nl 

VI 

19) 

9 

EFF1 

VI 

20) 

9 

VAX1 

VI 

21) 

9 

P2 

VI 

22) 

9 

P03 

VI 

23) 

9 

T2 

VI 

24) 

9 

T03 

VI 

25) 

9 

RH02 

VI 

26) 

9 

V2 

VI 

27) 

9 

BET2 

VI 

28) 

9 

ALF3 

VI 

29) 

9 

EFA2 

VI 

30) 

9 

P0W1 

VI 

31) 

9 

P3 

VI 

32) 

9 

T3 

VI 

33) 

9 

RH03 

VI 

34) 

9 

V3 

VI 

35) 

9 

CP2 

VI 

36) 

9 

KG2 

VI 

37) 

9 

S-R2 

VI 

38) 

9 

S-N2 

VI 

39) 

9 

EFF2 

VI 

40) 

9 

VAX2 

VI 

41) 

9 

P4 

VI 

42) 

9 

P05 

VI 

43) 

9 

T4 

VI 

44) 

9 

T05 

VI 

45) 

9 

RH04 

VI 

46) 

9 

V4 

VI 

47) 

9 

BET4 

VI 

48) 

9 

ALF5 

VI 

49) 

9 

EFA4 

VI 

50) 

9 

P0W2 

ITS  TRIAL  VALUE 
0.93567 
0.95267 
18.60127 
0.72132 
0.10240 
0.31627 
0.58133 
11.77661 
28.49974 
47.56911 
889.58203 
1040.01636 
4.53787 
916.76270 
2793.42432 
1.33821 
0.24222 
0. 14909 
0.82573 
254.05411 
20.39499 
21.48415 
871.95215 
883.72070 
3.31306 
256.41748 
1 .03000 
-0.13588 
0.00112 
0.41594 
16.53975 
820.62891 
2.85483 
588.54761 
2745.11304 
1.34623 
0.16165 
0.13204 
0.85735 
197.99052 
11.15817 
11.57824 
764.94360 
772.34985 
2.06615 
201.65057 
0.19062 
C. 19082 
0.00229 
0.29126 


0.0100 
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v< 

51) 

■ 

P5 

* 

9.14639 

V( 

52) 

* 

T5 

a 

722.31396 

v< 

531 

M 

RH05 

a 

1.79358 

vi 

541 

m 

V5 

a 

520.74292 

VI 

55) 

* 

CP3 

a 

2709.78857 

VI 

56) 

* 

KG3 

a 

1.35233 

VI 

571 

a 

S-R3 

a 

0.13131 

VI 

58) 

a 

S-N3 

r 

0.16165 

VI 

59) 

a 

EFF3 

a 

0.85791 

VI 

60) 

a 

VAX3 

at 

175.18066 

VI 

61) 

a 

P6 

■ 

6.46630 

VI 

62) 

a 

P07 

s 

6.67799 

VI 

63) 

m 

T  6 

■ 

678.06128 

VI 

64) 

a 

T07 

■ 

683.84497 

VI 

65) 

a 

RH06 

a 

1.35078 

VI 

66) 

a 

V6 

a 

177.04234 

VI 

67) 

a 

BE76 

a 

-0.14515 

VI 

68) 

a 

ALF7 

a 

-0.14515 

VI 

69) 

a 

EFA6 

a 

0.00398 

VI 

70) 

a 

POW3 

a 

0.22848 

VI 

71) 

a 

PEXO 

a 

6.62783 

VC 

72) 

a 

7EX0 

■ 

682.47412 

VI 

73) 

a 

RHOX 

V 

1.37557 

VI 

74) 

a 

VEXO 

■ 

85.39653 

VI 

75) 

a 

PTHT 

■ 

5.62490 

VI 

76) 

a 

77H7 

* 

579.20068 

VI 

77) 

a 

RHOT 

■ 

0.37836 

VI 

78) 

a 

V7H7 

■ 

746.10352 

VI 

79) 

a 

CPEX 

■ 

2659.82422 

VI 

80) 

a 

KGEX 

■ 

1.36134 

VI 

81) 

a 

ATH7 

a 

0.00337 
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■  KP.fc* 


VARIABLE  NUMBER  ANO  ITS  FINAL  VALUE 

THESE  VALUES  ARE  FOR  AN  ALTERNATOR  OUTPUT  OF  0.40  HM 


VARIABLE  VALUE  I  SI)  VALUE  I  BRIT ISH) 


VI  1) 

X 

POW 

X 

0.46784 

MW 

627.37939 

HP 

V(  2) 

= 

MFLO 

X 

0.47272 

KG/ SEC 

1.04217 

LBM/SEC 

VI  3) 

3 

MFLX 

3 

9.32930  KG/SEC-SQ.M 

1.91079 

LBM/SEC -SO. FT 

VI  4) 

x 

PHI 

S 

0.75766 

0.75766 

VI  5) 

s 

XNH3 

X 

0.08788 

0.08788 

VI  6) 

X 

XN2 

X 

0.31669 

0.31869 

VI  7) 

s 

XH2 

3 

0.59343 

0.59343 

VI  8) 

x 

MWT 

X 

11.62144 

KG/MOLE 

11.62144 

lbm/mole 

VI  9) 

X 

PI 

s 

13.92802 

8ARS 

202.00893 

PS!  A 

VI  10 ) 

2 

P01 

X 

23.30045 

BARS 

337.94434 

PSIA 

VIII) 

X 

T1 

X 

856.38135 

DEG  K 

1541.48560 

OEG  R 

VI 12 ) 

X 

T01 

X 

1004.46631 

OEG  K 

1808.03857 

OEG  R 

VI 13) 

X 

RH01 

X 

2.27331 

KG/CU.M 

0.14190 

LBM/CU.FT 

VI 14) 

X 

vi 

X 

908.05420 

M/SEC 

2979.18066 

FT/SEC 

VI 15) 

X 

CPI 

r 

2784.08521 

J/KG-OEG  K 

0.66541 

BTU/L8M-DEG 

R 

VI 16) 

x 

KGl 

3 

1.34584 

1.34584 

VI 17) 

8 

S-Rl 

3 

0.24222 

0.24222 

VI 18 ) 

X 

S-Nl 

X 

0.14909 

0.14909 

VI  19) 

X 

FFF 1 

3 

0.82710 

0.82710 

VI 20) 

X 

VAX1 

3 

251.64078 

M/SEC 

825.59302 

FT/SEC 

VI  21) 

X 

P2 

X 

9.88056 

8ARS 

143.30545 

PSIA 

VI  22) 

X 

P03 

X 

10.41240 

8ARS 

151.01915 

PSIA 

VI  23) 

e 

T2 

X 

837.55127 

OEG  K 

1507.59155 

OEG  R 

VI 24 ) 

s 

T  03 

X 

849.13599 

OEG  K 

1528.44409 

OEG  R 

VI  25) 

X 

RH02 

X 

1.64895 

KG /CU.M 

0.10293 

LBM/CU.FT 

VI 26) 

m 

V2 

X 

253.98177 

M/SEC 

833.27344 

FT/SEC 

VI  27) 

X 

BET  2 

3 

1.03000 

RAOIANS 

59.01897 

OEGREES 

V 1 28) 

X 

ALF3 

X 

*0. 1358B 

RADIANS 

-7.78581 

OEGREES 

VI  29) 

X 

EFA2 

X 

0.00113 

SQ.M 

0.01215 

SO. FT 

VI 30) 

X 

POWI 

X 

0.20443 

MW 

274.14478 

HP 

VI31) 

X 

P3 

X 

7.95697 

BARS 

115.40613 

PSIA 

V(  32) 

X 

T3 

X 

785.77026 

OEG  K 

1414.38574 

DEG  R 

VI33 ) 

X 

RH03 

3 

1.41543 

KG/CU.M 

0.08835 

LBM/CU.FT 

VI 34) 

3 

V3 

X 

589.02563 

M/SEC 

1932.49878 

FT/SEC 

V 135 ) 

3 

CP2 

3 

2737.67871 

J/KG-OEG  K 

0.65432 

BTU/LBM-OEG 

R 

VI 36) 

X 

KG2 

= 

1.35377 

1.35377 

V 137 ) 

X 

S-R2 

X 

0.16186 

0.16186 

VI 38 ) 

S 

S-N2 

X 

0.13204 

0.13204 

V 1 39 ) 

X 

EFF2 

X 

0.85725 

0.85725 

V 1 40) 

5 

VAX2 

X 

198.15137 

M/SEC 

650.10278 

FT/SEC 

VI41) 

X 

P4 

3 

5.29283 

BARS 

76.76602 

PSIA 

VI 42) 

3 

P05 

?: 

5.49942 

BARS 

79.76234 

PSIA 

VI43 ) 

3 

T  4 

3 

729.84863 

OEG  K 

1313.72681 

OEG  R 

VI 44) 

s 

T05 

X 

737.29297 

OEG  K 

1327.12671 

OEG  R 

VI45) 

3 

RH04 

3 

1.01366 

KG/CU.M 

0.06327 

LBM/CU.  FT 

VI46) 

X 

74 

X 

201.89430 

M/SEC 

662.38281 

FT/SEC 

V 147 ) 

X 

6FT4 

3 

0.19286 

RAO I  ANS 

11.05066 

OEGREES 

VI 48) 

X 

A’.F5 

X 

0.19286 

RAOIANS 

11.05066 

OEGREES 

VI49) 

X 

Ei"A4 

X 

0.00231 

SQ.M 

0.02486 

SO. FT 
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V(  50) 

X 

P0W2 

m 

0.14474 

MW 

194.10316 

HP 

V(  5  1 ) 

X 

P5 

8 

4.26603 

BARS 

61.87350 

PSIA 

V(  52) 

X 

T5 

X 

684.82446 

0E6  K 

1232.68335 

OEG  R 

VC53) 

s 

RH05 

9 

0.87073 

KG /CU.M 

0.05435 

LBM/CU.FT 

VC54) 

X 

V5 

9 

532.26147 

M/SEC 

1746.26465 

FT/SEC 

VC  55) 

X 

CP3 

X 

2699.73169 

J/KG-DEG  K 

0.64525 

BTU/LBM-OEG  R 

V(  56 ) 

X 

KG3 

X 

1.36054 

1.36054 

VC57) 

= 

S-R3 

9 

0.13637 

0.13637 

V(  58 ) 

X 

S-N3 

X 

0.16186 

0.16186 

VC  59) 

9 

EFF3 

X 

0.85595 

0.85595 

V 1  60) 

S 

VAX3 

X 

179.05553 

M/SEC 

587.45239 

FT/SEC 

V(  6 1 ) 

3 

P6 

9 

2.91084 

BARS 

42.21815 

PSIA 

V(62) 

* 

P07 

X 

3.01371 

BARS 

43.71017 

PSIA 

V(  63) 

8 

T6 

X 

638.33423 

OEG  K 

1149.00098 

DEG  R 

VC64) 

3 

T07 

9 

644.31201 

DEG  K 

1159.76099 

OEG  R 

V(  65) 

S 

RH06 

X 

0.63739 

KG /CU.M 

0.03979 

LBM/CU.FT 

VC66 ) 

X 

V6 

X 

179.66287 

M/SEC 

589.44507 

FT/SEC 

V)  67) 

9 

BET6 

X 

-0.08225 

RADIANS 

-4.71290 

DEGREES 

VI  68 ) 

9 

ALP  7 

X 

-0.08225 

RADIANS 

-4.71290 

OEGREES 

V)  69) 

9 

EFA6 

= 

0.00413 

SO.M 

0.04443 

SO. FT 

VC  70 ) 

9 

P0W3 

X 

0.11866 

MW 

159.13138 

HP 

VC  71 ) 

9 

PEXO 

9 

2.98756 

BARS 

43. 33090 

PSIA 

VC72) 

9 

TEXO 

X 

642.79736 

DEG  K 

1157.03467 

OEG  R 

VC  73 ) 

9 

RHOX 

n 

0.64965 

KG/CU.M 

0.04055 

LBM/CU.  FT 

VC74) 

9 

VEXO 

« 

B9. 72345 

M/SEC 

294.36816 

FT/SEC 

VC  75 ) 

3 

P7HT 

= 

2.52873 

BARS 

36.67613 

PSIA 

VC  76) 

9 

TTHT 

X 

544.07642 

OEG  K 

979.33691 

OEG  R 

VC77) 

X 

RHOT 

s 

0.18210 

KG/CU.M 

0.01137 

LBM/CU. FT 

VC78) 

9 

VTHT 

8 

729.84009 

M/SEC 

rr\ 

(A 

00 

CO 

• 

* 

O' 

<<\ 

rvj 

FT/SEC 

VC79J 

X 

CPEX 

X 

2657.07178 

J/KG-OEG  K 

0.63506 

BTU/LBM-DEG  R 

VC  80) 

9 

KGEX 

X 

1.36846 

1.36846 

■VC  81) 

9 

ATHT 

9 

0.00356 

SO.M 

0.03829 

SO.  FT 
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APPENDIX  C 

The  particular  input  formats  employed  in  this  program  are  explained 
at  the  beginning  of  APPENDIX  B.  Example  data  cards  for  a  2.0-mw,  three- 
stage  turbine  system  follow: 


Data 

Card 

Number 


Colunn  Number 


1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

61-70 

71-80 

2.4 

2.4 

47.4 

0.66 

0.13 

0.31 

0.56 

12.0 

74.5 

124. 

964. 

1120. 

11.2 

939. 

2824. 

1.32 

0.24 

0.15 

0.82 

260.0 

54.4 

57.2 

950. 

962. 

8.29 

262. 

1.03 

-0.13 

0.001 

1.07 

44.8 

899. 

7.2 

588.6 

2774. 

1.33 

0.16 

0.13 

0.86 

198. 

31.1 

32.2 

844.6 

851.9 

5.33 

201.7 

0.19 

0.19 

0.002 

0.735 

26.3 

806.6 

4.7 

499.5 

2750. 

1.33 

0.12 

0.16 

0.86 

168. 

19.8 

20.4 

767. 

772. 

3.7 

74.1 

-0.26 

-0.26 

0.004 

0.53 

20.2 

771. 

3.8 

78.0 

17.2 

658. 

1.00 

782. 

2682. 

1.35 

0.003 

eBM 

5-8 

9-12 

13-16 

17-20 

21-24. . 

• 

..77-80 

